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Résumé en Anglais 


This work is devoted to the study of the existence, uniqueness, periodicity 
and positivity of solutions for two classes of delay second order differential 
equations. The first class is of neutral differential equations with time varying 
delay while the second one is of iterative differential equations with time and 
state dependent delays. The method used here is one of the most efficient 
techniques for studying these types of equations since it combines some useful 
properties of Green’s functions together with fixed point theorems to esta- 
blish sufficient conditions for proving the desired results. The idea of this 
technique is based on the converting of the considered equation into an inte- 
gral one whose solutions are recourse to a Krasnoselskit’s or Schauder’s fixed 
point theorem. In addition, by applying the principle of contraction mapping 
on the integral equation and under certain specified constraints we proved 
the existence and uniqueness of the periodic solution of the second class of 


equations. 


Keywords : Contracting principle, existence, Green’s function, iterative dif- 
ferential equation, Krasnoselskii’s fixed point theorem, neutral differential 


equation, Schauder’s fixed point theorem, uniqueness. 


Résumé en Francais 


Ce travail a été dévoué a l’étude de l’existence, l’unicité, la périodicité et 
la positivité des solutions de deux classes d’équations différentielles du second 
ordre à retard variable. La première classe est d'équations différentielles de 
type neutre avec un retard variant dans le temps tandis que la deuxième est 
d'équations différentielles itératives avec des retards dépendants du temps et 
de l’état. La méthode utilisée ici est l’une des techniques les plus efficaces pour 
étudier ces types d'équations, car elle combine certaines propriétés utiles des 
fonctions de Green avec des théorèmes de point fixe afin d’établir des condi- 
tions suffisantes pour prouver les résultats souhaités. L’idée de cette technique 
est basée sur l’inversion de l’équation considérée en une équation intégrale 
dont la solution fait recoure au théorème de point fixe de Krasnoselskii ou 
celui de Schauder. Aussi, moyennant le principe de l’application contractante 
de Banach et sous certaines contraintes déterminées, on a montré l’existence 


et l’unicité de la solution de la deuxième classe d'équations. 


Mots clés : Principe de contraction, existence, fonction de Green, équa- 
tion différentielle itérative, théorème de point fixe de Krasnoselskii, équation 


différentielle de type neutre, théorème de point fixe de Schauder, unicité. 
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Po CT pographie et Acronymes 


Typographie 
T désigne le retard dans l’ensemble de cette thèse et est, soit une fonction, 
soit une constante suivant les modèles et les exemples. 
u désigne l’argument dévié dans l’équation à arguments déviés. 
T désigne la période. 
Les titres des références, les noms des savants, les termes techniques ou les 
mots que nous voulons souligner sont indiqués en italique. 
Les notations entre guillemets correspondent à l’emploi de termes techniques 


que nous voulons faire ressortir. 


Acronymes 


Abréviation Signification 


CSH cellule souche hématopoïétique. 
dde23 solveur de Matlab, permet de résoudre les EDR à retard constant. 
ddensd solveur de Matlab, permet de résoudre des équations différentielles, 


de type neutre où le retard est indépendant de l’état. 
ddesd solveur de Matlab, permet la résolution des équations différentielles, 


à retard variable. 


Acronymes et Typographie 


ddefile 
ddefun 
DELSOL 


EDF 
EDN 
EDO 
EDP 
EDPR 
EDR 
DELSOL 


Fortran 


history 


lags 


Maple 


Matlab 


P-P 
tspan 


RADARS 


RETARD 


fonction au second membre de l’équation différentielle à retard. 
fonction au second membre de l’équation différentielle à retard. 
solveur de Fortran, permet de résoudre les équations différentielles, 
à retard. 

équation différentielle fonctionnelle. 

équation différentielle de type neutre. 

équation différentielle ordinaire. 

équations aux dérivées partielles. 

équations aux dérivées partielles à retard. 

équation différentielle à retard. 

solveur de Fortran, permet de résoudre les équations différentielles, 
à retard. 

langage de programmation compilé utilisé principalement pour, 

le calcul scientifique. 

fonction historique. 

vecteur colonne représente les retards. 

langage interprété de programmation utilisé principalement pour, 
le calcul formel. 

langage de programmation interprété utilisé principalement pour, 
le calcul matriciel et numérique interactif 

Presque partout. 

intervalle d'intégration. 

solveur de Fortran, permet de résoudre les équations différentielles, 
à retard. 

solveur de Fortran, permet de résoudre les équations différentielles, 


à retard. 


Pe Notations 


R”)) 


Ensembles et nombres 


ensemble des nombres réels. 

ensemble des nombres réels positifs ou nuls. 
espace euclidien de dimension n. 

ensemble des entiers naturels. 


ensemble des entiers naturels non nuls. 


intervalle fermé de R d’extrémités a et b. 


intervalle semi-ouvert de R d’extrémités a et b. 


intervalle ouvert de R d’extrémités a et b. 


ee 


espace des fonctions continues de |a, b] dans 
espace de Sobolev d’ordre 1. 

espace de Lebesgue des fonctions sommables, 

sur (0, +00). 

ensemble des fonctions continues et T— périodiques. 
ensemble des nombres complexes. 

ensemble des fonctions continues et T'—périodiques. 
ensemble des fonctions continues et T— périodiques. 


telles qu’il existe K >Oet L>0:K<v< L. 


Xi 


Notations 


Pr (L, M) 


| 

II Flloo 
max 
min 
sup 
expr 


G (t,x) 


ensemble des fonctions x de Pr de norme ||zx|| < L, 


telles que |z (t2) — x (t1)| < M |te — til, Vti, te ER. 


r= T est la dérivée de la variable x par rapport au temps et, 
c’est la dérivée à droite dans les résultats d’existence et d’unicité 
dans la deuxième section du deuxième chapitre. 

variable temporelle. 

boule fermée de centre a et de rayon r. 


espace de toute les suites réelles avec lim Ç; = 0. 


1—00 


K=R ou C. 


Fonctions 


valeur absolue. 

norme de la convergence uniforme et se définit par sup |f (x)|. 
fonction maximum. 

fonction minimum. 

fonction supremum. 

fonction exponentielle de x. 


fonction de Green. 
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Introduction générale 


D” la nature, plusieurs phénomènes sont gouvernés par une classe 
d'équations différentielles à retard. Comme leur non l'indique, ces 
équations notées en abrégé EDR modélisent généralement l’évolution des va- 
riables dépendant non seulement des valeurs actuelles mais dépendent aussi 
irréductiblement des valeurs prises dans le passé autrement dit elles tiennent 
compte de l’effet du passé dans la prédiction du futur c’est pour ça qu’elles 
sont souvent appelées "les équations à mémoire", "les équations héréditaires", 
"équations à post-effet", ou encore "équations à argument dévié". 

Le retard est généralement une constante positive, une variable dépen- 
dante continiment du temps ou de l’état ou distribué et se traduit comme 
un temps nécessaire pour que le système réponde à une certaine évolution, 
ou parce qu’un certain seuil doit être atteint avant que le système ne soit 
activé. 

Sa signification diffère d’un modèle à un autre, il peut représenter la durée 
du cycle cellulaire dans la division cellulaire, le temps de gestation ou de 
développement dans la dynamique des populations, la période d’incubation 
d’une maladie contagieuse en épidémiologie etc. 


L'introduction d’un tel retard a fait cauchemarder plusieurs chercheurs 
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pendant un siécle car sa présence peut induire des oscillations indésirables 
ou peut causer l’instabilité du système et rend la théorie difficile et conduit 
à faire appel à un arsenal mathématique très vaste. Malgré ça, il a fait rêver 
d’autres aussi comme Picard qui a souligné lors de la conférence interna- 
tionale des mathématiciens à Rome en 1908, l’importance de modéliser les 
systèmes physiques par des équations retardées et il a affirmé : 

Les équations différentielles de la mécanique classique sont telles qu'il en 
résulte que le mouvement est déterminé par la simple connaissance des posi- 
tions et des vitesses, c’est-à-dire par l’état à un instant donné et à l’instant 
infiniment voisin. 

Les états antérieurs n’y intervenant pas, Vhérédité y est un vain mot. 
L'application de ces équations où le passé ne se distingue pas de l’avenir, où 
les mouvements sont de nature réversible, sont donc inapplicables aux êtres 
vivants. 

Nous pouvons rêver d'équations fonctionnelles plus compliquées que les 
équations classiques parce qu’elles renfermeront en outre des intégrales prises 
entre un temps passé très éloigné et le temps actuel, qui apporteront la part 
de l’héréditä] 

D'autre études ont montré également que l’introduction volontaire du 
retard peut améliorer la stabilisation des systèmes : amortissement et stabili- 
sation [I], résonateurs retardés [45], rejet de perturbation [46], [89], contrôle 
de cycle limite non linéaire [4]. 

La littérature concernant les équations fonctionnelles à retard est abondante, 
on ne peut pas citer le grand nombre d'ouvrages dédiés à cette thématique 
mais le lecteur intéressé pourra se référer aux livres de Hale-Cruz |40], Hale 
[41], Hale- Verduyn Lunel [42], Kolmanovwskii-Myshkis A9] et Smith [BI]. 


'J.K. Hale, History of delay equations, dans [9], 1 — 28. 
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Plusieurs approches ont été utilisées pour traiter ce type d’équations 
comme la méthode des étapes, la méthode de Runge Kutta, les transfor- 
mations de Laplace et de Schröder et la théorie de point fixe qui est un 
outil précieux pour la résolution des équations différentielles non linéaires en 
transformant le probléme d’existence en un probléme de point fixe. 

Le fil rouge de cette thése est d’appliquer cette derniére technique pour 
étudier certaines classes d’équations différentielles non linéaires à retard lar- 
gement rencontrées dans la modélisation de plusieurs phénomènes dans des 
domaines très variés comme la physique, la chimie, la biologie, l’écologie, 
etc. Plus précisément, nous nous intéresserons à appliquer une technique qui 
combine un des théorèmes de point fixe de Banach, de Schauder ou de Kras- 
noselskii avec quelques propriétés utiles des fonctions de Green de certaines 
équations différentielles linéaires pour arriver aux résultats désirés. 

Donc, notre travail de thèse s'intéresse aux problématiques suivantes : 
l'existence, l’unicité, la périodicité et la positivité des solutions de deux 
classes d'équations différentielles à retards variables, la première est de type 
neutre avec des retards variant dans le temps et la deuxième est de type 
itérative avec des retards dépendants de l’état. 

On a structuré ce manuscrit en quatre chapitres et une annexe sur la syn- 
taxe des principaux solveurs de Matlab et Maple qui permettent la résolution 
numérique des équations différentielles à retard. 

Dans le premier chapitre, nous avons jugé utile de rappeler brièvement 
quelques notions générales et certains résultats fondamentaux indispensables 
pour le reste de la thèse comme la fonction de Green unidimensionnelle et 
le théorème d’Ascoli-Arzela où la partie principale du chapitre est dédiée 
à la présentation des théorèmes de point fixe, le théorème de point fixe de 


Banach, de Schauder et celui de Krasnodelskii. 
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Le deuxième chapitre est un état de l’art sur les équations différentielles à 
retard. Après un bref historique et une classification de ce type d’équations, 
nous passons à un rappel de quelques notions de base et certains résultats de 
la théorie d’existence et d’unicité. Ensuite, nous nous intéresserons plus par- 
ticulièrement à présenter quelques méthodes de résolution comme la méthode 
des étapes, la méthode de Runge-Kutta et la transformation de Laplace. Nous 
proposerons pour conclure des modèles à retard modélisant cinq phénomènes 
dans l’écologie, l’électrodynamique classique, la biologie et l’électronique. 

Enfin, l'essentiel de cette thèse est présenté dans le troisième et le qua- 
trième chapitres qui correspondent à des articles publiés [15] [16] ; seules les 
références ont été regroupées dans une bibliographie commune. La méthode 
utilisée dans ces deux chapitres est très efficace car elle combine certaines 
propriétés utiles des fonctions de Green et la technique du point fixe. Dans 
le troisième chapitre, on va utiliser le théorème de Krasnoselskit pour établir 
l'existence des solutions périodiques positives de l’équation différentielle non 


linéaire du second ordre de type neutre suivante : 


ult) + rll) Lat) + a(or() = a (t, z(t — r(t))) 


+f(t, hı (x(t), (w(t — T(t)))), 


où p et q sont des fonctions positives continues. La fonction Q : R x R — R 


est différentiable, h1, gı : R — R et f : R x R x R — R sont des fonctions 
continues par rapport à leurs arguments. 

Le travail exposé dans ce chapitre est réalisé en collaboration avec Pr 
Ahcene. Djoudi. et Dr Abdelouaheb. Ardjouni., et est publié dans le journal 
Nonlinear Dynamics and Systems Theory [15]. 


Le quatrième et dernier chapitre est dédié à l’étude de l’équation diffé- 
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rentielle itérative du second ordre suivante : 


Sut + p (t) © x(t) talt) e(t) = “9 (t,2(t), xl), ..., 2(2)) 


+h (t, e(t), sP (t), c(t) , 


où z(t) = t, cA) = a(t), A(t) = z (x(t), s(t) = cl (x(t)). Les 
résultats exposés dans ce chapitre est une publication dans le journal Acta 
Mathematica Universitatis Comenianae [16]. 

On termine cette thèse par une conclusion générale sans doute non ex- 


haustive ainsi que les perspectives de nos recherches futures. 
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A cause de leur grande importance dans notre étude, ce chapitre intro- 
ductif vise à présenter quelques notions de base nécessaires pour le 
développement ultérieur de cette thèse et on insistera en particulier sur les 
fonctions de Green unidimensionnelles, au théorème d’ Ascoli-Arzela et aux 


théorèmes de point fixe de Banach, de Schauder et celui de Krasnoselskii. 


Chapitre 1. Notions Préliminaires 


1.1 Fonctions de Green uni-dimensionnelles 


La fonction de Green] qui est une solution élémentaire d’une équation 
linéaire que ce soit différentielle ou aux dérivées partielles, est un outil très 
puissant pour la résolution de ces équations et son intérêt principal réside 
dans la représentation intégrale de la solution. 

On va aborder cette notion que très légèrement ici, juste pour rappeler 
les conditions de son existence. 
ième 


Soit l’équation différentielle du n ordre suivante : 


(n—1) 


ao (t)y™ (t) tar (ty (+. tan (t) y(t) = A(t), Vte[e,d], (11) 


où les fonctions ao, 1, ..., Gn, h sont continues sur |c, d] , ao (t) 4 0 sur [c,d] 


avec les conditions aux limites suivantes : 


any (c) + ay (ec) +... FA Vy (6) + Bay (d) +... + BEPy (d) = 0. 
(1.2) 


L’équation homogène de (1.1) est : 


(n—1) 


ao (t)y™ (t) + a lty (+. tan (t)y(t) =0, yte led. (13) 


‘Cette fonction doit son nom au physicien britannique George Green (1793 — 1841). 
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On suppose que le rang du tableau A ci-dessous est n : 


ay a) se Di m ge» 
anle aP u Ba … BET 
G Oe ai P ae Y 


Théorème 1.1 [60] 

a) Si l'équation homogène de avec les conditions frontières nad- 
met que la solution triviale, alors il existe pour ce probléme une fonction de 
Green unique G (t, s) on suppose ao (t) £ 0, Vt € [c,d]. 

b) on suppose réalisées les conditions ci-dessus, alors l’équation non homo- 
gène admet une solution unique vérifiant les conditions frontières 
donnée par : 


y (t) -f G (t,s) g(s)ds. (1.4) 


Supposons que le problème aux limites (1.3)-(1.2) n’admet que la solution 
triviale y (t) = 0. 


Définition 1.1 [52] On appelle "fonction de Green" ou "fonction d'influence" 
de l’équation (1.3) avec les conditions frontières (1.2), la fonction G (t, s) 
construite pour tout point s € |c,d[, et jouissant des quatre propriétés sui- 
vantes : 

1/ G (t, s) est continue et possède des dérivées continues par rapport à t jus- 
qu’à l’ordre (n — 2) inclus pour t € [c, d]. 


2/ Sa (n — j dérivée par rapport à t présente au point t = s une discon- 


1 
tinuité de première espèce, le saut ayant la valeur (s)’ i.e. 
ag\s 
n—1 n—1 
o ig o a j= 1 
on | or! ao (s)’ 
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ou 

c (st, s) et “c (57, s8) 
désignent respectivement les limites à droite et à gauche lorsque t tend vers 
s de la fonction qui à t fait correspondre (t, s). 
3/ Pour tout s fixé dans |c, d|, cette fonction vérifie la condition frontière en 
c et la condition frontière en d. 
4/ Pour tout s fixé dans|c, d| , la fonction qui à t fait correspondre G (t, s) est 


solution de l'équation homogène dans chacun des intervalles [c, s) et (s, d]. 


1.2 Théorème d’Ascoli-Arzela 


Ci-dessous on rappelle le théorème d’ Ascoli-Arzela qui est un outil classique 
et puissant pour montrer qu’une partie de l’espace des fonctions continues 


sur un compact est relativement compacte. 


Définition 1.2 Soit M un sous ensemble de E = C ([a,b], R), M est dit 


équicontinu si et seulement si 


Ve > 0, Vz € [a,b], 3ô > 0, Vy € [a,b], 
[lz — yll <= Nf € M, |f(x) — f(y)| < el. 
Exemple 1.1 


Si k est un nombre réel positif, alors l’ensemble des fonctions lipschitziennes 


de E dans R de rapport k est équicontinu. En effet, pour e fixé, il suffit de 


prendre 

8 
M ? 

dans la définition 1.2. 
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Théorème 1.2 [90] Soit C([a,b],R), —oo < a < b < +, l’espace des 
fonctions continues définies sur le compact |a, b] et à valeurs réelles muni de 
la norme 


= t)|. 
lu = max lue) 


Une partie F de C ([a,b],R) est relativement compacte dans C (la, b] ,R) si et 


seulement si elle est uniformément bornéd?] et équicontinue. 


Remarque 1.1 Le théorème d’ Ascoli-Arzela ne permet de caractériser que 


les ensembles relativement compacts de C (E4, K), (avec E compact et mé- 


trique, K = R ou C), et non les ensembles relativement compacts de n’im- 


porte quel espace de Banach. 


Exemple 1.2 Soient kı et kz deux réels strictement positifs. Le sous-ensemble 


F des fonctions réelles continues sur [a,b], dérivables sur |a, b| qui vérifient 


|f Æ| < kı et sup |f (t)| < kə, 


pour tout t € [a,b], est relativement compact dans C (|a, b], R). 
En effet, pour tout f € F, le théorème des accroissements finis, prouve que 


pour tout to, t € [a,b] il existe c € to, t| tel que 


IF) — f (to)| = |F (| |t — tol - 
Donc |f (t) — f (to)| < kalt — tol. Fixons to € [a,b]. Soit € > 0 et n = a 
2 
alors 


vt € [a,b], |t— tol < n = |f &) — F (to)| Se. 


2F est dite "uniformément bornée" si et seulement si 


Va € [a,b], {f (x), fE F} est borné. 
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Ce qui est exactement l’équicontinuité de F en ty. Comme nous pouvons 
prendre pour to n’importe quel point de [a,b], on en déduit que F est équi- 


continu. 


On a |f (t)| < kı pour tout f € F ce qui implique que || f||,, < kı et partant 
Yf EF, feB'(0,k), 


i.e., 


Fe B'(0, kı), 


d’où la bornitude de F. 
Finalement, Comme F est borné et équicontinu, alors le théorème d’ Ascoli- 


Arzela assure que F est relativement compact. 


1.3 Théorémes de point fixe 


Comme l’attestent les très nombreux travaux paraissant aujourd’hui au ni- 
veau international, les théorèmes de point fixe sont des outils précieux et 
très intéressants en mathématiques, surtout pour la résolution des équations 
différentielles non linéaires. Pour résoudre un problème par la technique du 
point fixe, nous avons besoin d’une application appropriée, d’un ensemble 
convenable apte pour contenir les solutions du problème et d’un théorème de 
point fixe qui donne certaines conditions sous lesquelles cette application ad- 
met au moins un point fixe. On va voir maintenant trois théorèmes de point 
fixe, le théorème de point fixe de Banach qui donne un critère général dans 
les espaces métriques complets, celui de Schauder qui est plus topologique et 
affirme qu’une application continue sur un convexe compact admet au moins 


un point fixe et finalement, le théorème hybride de Krasnoselskii. 
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Définition 1.3 Soit (F, ||.||,,) un espace de Banach, T : E — E une ap- 
plication. On appelle point fixe de T tout point x € E tel que T(x) = x. Ce 


qui est équivalent à dire que l’équation T(x) — x = 0 possède une solution. 


Définition 1.4 Soit (E, ||.||,,) et (F, ||.) deux espaces de Banach et 
T : E — F une application. On dit que T est lipschitzienne de rapport 
k > 0 sur E si 


IT -Tyle < k lle —ylle- (1.5) 


En particulier, si 1 > k > 0, T est dite contraction ou application contractante 


de rapport k. 


1.3.1 Théorème de point fixe de Banach 


En 1922, le mathématicien polonais Stefan Banach a prouvé son célèbre 
théorème (connu aussi sous le nom de théorème de l'application contractante 
ou théorème de Banach- Picard) qui garantit l’existence et l’unicité d’un 
point fixe d’une application contractante d’un espace métrique complet dans 
lui-même. En outre, il est basé sur un processus itératif assurant que ce point 
fixe peut être obtenu comme limite d’une suite itérée et qu’il est possible 


d'estimer la précision avec laquelle cette limite est atteinte. 


Théorème 1.3 [78] Soit (E, ||-||g) un espace de Banach et T : E — E 
une application contractante de constante k € [0,1[. Alors il existe un point 


unique x € E tel que T(x) = x. 


Remarque 1.2 Si À est une application lipschitzienne (pas nécessairement 


une contraction) mais l’une de ces itérées A? est une contraction, alors A a 


3 Charles Emile Picard fut le premier à utiliser le théorème de point fixe de Banach 


dans une méthode d’approximations successives. 
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un seul point fixe. 


En effet, soit x l’unique point fixe de A? on a 
AP (A(x)) =A (AP (x)) = À (x), 


ce qui convient à dire que A (x) est aussi un point fixe de A? et grâce à 
Punicité À (x) = x. 
Ce résultat est valable pour tous les types de contractions qui assurent l’uni- 


cité du point fixe. 


1.3.2 Théorème de point fixe de Schauder 


Le théorème de point fixe de Schauder élaboré en 1930, assure l’exis- 
tence d’au moins un point fixe pour une application continue sur un convexe 


compact dans un espace de Banach. 


Théorème 1.4 [78] Soit M un sous ensemble convexe, fermé, borné et non 
vide d’un espace de Banach E et T : M — E une application compacte. 


Alors T possède un point fire. 


Remarque 1.3 Si M est compact et convexe, il suffit que T soit continue 


pour avoir un point fixe pour T. 


1.3.3 Théorème de point fixe de Krasnoselskii 


En 1954, Kranoselskii combina le théorème de point fixe géométrique de 
Banach et le théorème de point fixe topologique de Schauder en un théorème 
hybride qui affirme que dans un convexe compact, toute application qui s'écrit 
sous forme d’une somme de deux applications dont l’une est une contraction 


et l’autre est compacte admet au moins un point fixe. 
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Principe du théorème de point fixe de Krasnoselskii [17] 
Selon Smart [78], Krasnoselskii étudia un article de Schauder [74] et formula 
le principe suivant ; 

L’inversion d’un opérateur différentiel perturbé peut engendrer une somme 
de deux opérateurs, une contraction et un opérateur compact. 

Pour mieux comprendre cette observation de Krasnoselskii, on considère 


l'équation différentielle suivante : 
a(t) = —a(t)a(t) — g(t, x(t), (1.6) 


où a(t +T) = a(t) et g(t +T, x) = g(t, x) pour un certain T > 0. 


On peut transformer cette équation sous une autre forme en écrivant, 


x’ (texp ( | ads) = ~a(t)x(t)exp ( / ads) —g(t, «(t))exp ( f ads) 


(1.7) 
Par conséquent, 


(oDe f ass) does ( I ads) | (1.8) 


Une intégration de t — T à t donne 


[. (2twyern f aljas) du = -f ousa) (er f aas) du. 


(1.9) 


Ainsi, 


z(t) = z(t — T)ezp (- fos) z / | gluatu) exp (- i ads) du. 


(1.10) 
a (- | “(sds eae 
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et si (E, ||.||,,) est l’espace de Banach des fonctions y : R — R continues et 


T—périodiques, alors l’équation (1.10) peut se mettre sous la forme 


y(t) = (Be) (t) + (Ay) t) C= Pelt), (1-11) 


où B est une contraction de constante a < 1 et A est une application com- 
pacte. Cet exemple montre bien la naissance de l’application Py := By + Ap 
qui s’identifie à une somme d’une contraction et une application compacte. 

La recherche d’une solution pour exige donc un théorème adéquat qui 
s’applique à cet opérateur hybride P et qui peut conclure l’existence d’un 
point fixe qui sera à son tour, solution de l’équation initiale (1.6). Krasno- 
selskii trouva la solution en combinant les deux théorèmes de Banach et celui 
de Schauder en un seul théoréme hybride mais puissant qui porte son nom. 


Il établit le résultat suivant : 


Théorème 1.5 [78] Soit (E, ||.||,,) un espace de Banach et soit M une partie 
convexe, fermée et non vide de E. On suppose que À, B : M — E sont deux 
applications satisfaisant : 

(i) Ax + By € M, Yz, y E€ M; 

(ii) A est continue et A est compacte ; 

(iii) B est une contraction de constante a < 1. 


Alors, il existe x* € M tel que 
Ax* + Ba* = x*. 


Preuve. D’une part, la 3°"° condition donne 


IO —B)(x) —(1—B)(y)l = I(x- y) - (Br - By)| 
< x — yll + Bx — Byll 
< |e —yll+ alla -yl 
< (1+a)lx- yl, (1.12) 
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et d’autre part 


I —B)(x)-(—B){(y)} = I(@—y) — (Ba — By) | 
= | (@—y)|| — Ba — Byll 
> |æ-yll -alx -— yll 


> (1-=a)ļlz- yl. (1.13) 


De et (1.13), il résulte que 
(1 — a) ||x — yl] < A=B (z) -( - 8B) (y) || <(1+a)lx- y]. (114) 


Cette double inégalité montre que (J — B) : M — (I — B) M est continue 
et bijective. Donc, (Z — B) existe et elle est continue. 

Posons U := (I — B) A, il est clair que U est une application compacte, 
puisque U est une composition d’une application continue et une application 
compacte. 


En vertu du théorème de Schauder, U admet un point fixe, i.e. 
3 x € M tel que (1—B) Ar = x. (1.15) 
Ceci équivaut à dire 
Az + Br = 7, (1.16) 


et la preuve est achevée. m 


Remarque 1.4 Si A = 0, le théorème se résume au théorème de Banach et 


si B = 0 alors le théorème n’est autre que le théorème de Schauder. 


Remarque 1.5 [6I] Dans une des versions généralisées de ce théorème, la 
condition (iii) est remplacée par la condition suivante : Il existe une fonction 
continue à droite 


® : [0,+oo) — [0,+00), 
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telle que ¢(r) < ¢(s) pour 0 <r < s, ġ(r)<r sir > 0, et 
[Bx — By|| < ¢ (|x — yll) pour tout z, y € M. 


Remarque 1.6 Dans le cas où a = 1 dans la condition (iii), il n’y a 


pas de point fixe. En effet, voici un contre-exemple. 


Soit l’espace de Banach Co (R) de toutes les suites réelles avec lim ¢; = 0; 


1—00 


1:86, 


Co(R) = 16 = (Chien, GER: lim ¢, =o}, 


i—> oo 


muni de la norme 


Ita] = max {|¢;| :i = 1,2, nr ¢ = (Chien € Co ( R). 


On considère l'application A définie par 


A: B'(0,1) — OR) 
C= (Gien > AC=n= (i) scn ; 


où B'(0,1) désigne la boule unité fermée de Co (R) et 
1 = 2 (1 + 1101), 7 = A - ial) pour à = 2,3,... 
Puisque |7,| < 1 et |7,;| < |¢;| < 1, alors, on obtient 
A: B'(0,1) — B'(0,1). 


En outre, si ¢ = (C);ew et C= (ci) sont deux éléments distincts de 
iEN* 
B'(0,1), on aura 


AÇ- AG =n-1 = (n, - De 


où 


M =m = 27 (1 + (NCI) = 27 + EN) = 2746 = EN), 
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et 


nN; = (1-2 16, (1-2 YG, = (1-27) (ci = ea) pour 2 = 2,3,... 


Alors 


LAC — AC = max|n, —7,|:i =1,2,... 
= max {2 (6l = IC), max {(1- 2°) [Ga — Gl:5=2,8,.}} 
< ||¢-¢l. 


S’il existe Ç € B’(0,1) tel que AÇ = Ç alors, €, = 271(1 + |¢||) > 0. De plus 


pour tout 2 > 2, 


ICs | = a — D = a =- 20 = je 
i—2 i—2 
D i 
k=0 k=0 


i+1 3 
i = > (Gil > 1 Gil, 
j=3 
ce qui contredit le fait que Ç; tend vers zéro quand 7 tend vers linfini et 


prouve que A n’a pas de point fixe. 


1.4 Récapitulatif 


Dans ce chapitre, nous avons abordé brièvement quelques rappels et no- 
tions de base en vue de les utiliser dans les chapitres suivants comme la 
fonction de Green et le théorème d’Ascoli-Arzela qui donne les conditions 
suffisantes pour qu’une partie d’un espace de Banach soit relativement com- 
pacte sans oublier de présenter les trois théorémes captivants de point fixe, le 
théorème de Banach, de Schauder et celui de Krasnoselskii qu'ils ont initié de 
nombreuses études et ont été étendu dans plusieurs directions en modifiant 


certaines de leurs hypothèses. 
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es objectifs de ce chapitre sont multiples. Après un aperçu historique 
L. les EDR, la deuxième section sera consacré entièrement à la classi- 
fication de ces équations retardées selon les retards cités dans la littérature 
avant de nous intéresser à résoudre quelques EDR par trois méthodes : la mé- 
thode des étapes et de Runge-Kutta et la transformation de Laplace. Enfin, 


nous avons trouvé judicieux de présenter certains modèles retardées. 


Chapitre 2. Equations différentielles à retard 


2.1 Contexte historique 


Les équations fonctionnelles à retard constituent un champ d’étude très 
important pour modéliser des phénomènes d’hérédité rencontrés en physique, 
biologie, chimie, économie, écologie, etc. Malgré que dans la plupart des mo- 
dèles, le retard est estimé non signicatif et ignoré pour simplifier l’étude, il a 
été prouvé que dans de nombreux cas, le retard joue un rôle dominant dans 
plusieurs domaines et que les modèles avec retard fournissent des résultats 
plus précis et réalistes que leurs homologues sans retard. 

À notre connaissance l'apparition de ces équations remonte au 18°” 
siècle, elle est due à J. Bernoulli, L. Euler, J.L. Lagrange, P. Laplace, S. 
Poisson et d’autres. Par exemple, en 1728 dans ses expériences sur la corde 
vibrante et en partant d’une équation aux dérivées partielles de type hyper- 


bolique, Bernoulli [TA] a trouvé l'équation à retard suivante : 
y =y(t—1), 


mais malheureusement il décida la tenir pour fausse et il est dit qu’il y a eu 
plusieurs erreurs à déduire l'équation. 

Ce type d'équations a resté terra incognita jusqu’au début du 20 siècle 
et les travaux pionniers qui établissent le début de la théorie ont été dans 
la géométrie et la théorie des nombres et les premiers papiers traitant les 
équations fonctionnelles retardées linéaires sont dus à Polossuchin (1910) 
[67] et Schmidt (1911) [75]. 

Dans l’âge d’or de l'écologie théorique, l’étude de ce type d’équations 
connut un essor considérable avec les séries de travaux de V. Volterra 
sur les modèles prédateur-proie et les modèles de viscoélasticité. Il a 
utilisé la méthode d'énergie pour étudier une classe générale d'équations à 


retard non linéaires et il a écrit dans son œuvre majeure sur le rôle des 
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effets héréditaires sur les modèles de la dynamique de plusieurs espèces en 
interaction. 

Il y a eu d’intensives recherches sur le sujet depuis 1940 (surtout en lex- 
Union soviétique). La régulation sur base de modèles linéaires et stationnaires 
avec retard fut abordée en 1941 par Y. Zypkin. En 1949, A. D. Myshkis 
a posé les bases de la théorie moderne des EDR. En particulier, il fut le 
premier à formuler l’énoncé du problème de Cauchy pour des équations à 
retard arbitraire. 

Les années cinquante, ont vu une explosion de la théorie qui a été large- 
ment développée et les EDR fait partie du vocabulaire des chercheurs tra- 
vaillant sur la viscoélasticité, les problèmes mécaniques, les réacteurs nu- 
cléaires, le flux de chaleur, les réseaux de neurones, la combustion, l’interac- 
tion des espèces, les modèles microbiologiques, épidémiologiques ou physio- 
logiques, ainsi que beaucoup d’autres (voir [49]). 

La littérature concernant les EDR dans cette période est abondante et de 
nombreux travaux ont établi des résultats génériques dont on cite les résultats 
de Krasovskii [53] [54], qui a étendu la deuxième méthode de Liapunov aux 
équations retardées, sans oublier aussi les travaux de Myshkis [62], Krasovskii 
54], Bellman et Cooke [12], Halanay [39]. 

Les années suivantes ont donné naissance à un grand nombre de travaux 
dans cette direction, et surtout ceux qui concerne l’analyse de la stabilité 
des équations différentielles, avec un argument retardé comme El’sgol’ts et 
Norkin [33], Hale [4I], Hale et Lunel [42], Diekmann, Van Gils Lunel et 
Walther [22]. 
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2.2 Classification 


Le but de cette section est d’esquisser un tour d’horizon des principaux 
types d’équations différentielles 4 retard. Généralement, la pratique de la mo- 
délisation prouve que seules les équations retardées ou de type neutre sont uti- 
lisées pour représenter des phénomènes réelles mais il existe un nombre consi- 
dérable de phénomènes modélisées par des équations avancées, par exemple 


l'interprétation physique de l'équation différentielle avancée suivante : 
e(t)=2"(¢+1), 


est celle d’une force liée par exemple à l’élongation d’un ressort et qui n’agirait 
qu'après un certain temps. 

On a entamé cette section par la présentation de la classe mère dite des 
équations différentielles à argument dévié avant d’exposer les différents types 


d'équations différentielles à retard. 


2.2.1 Equations différentielles à argument dévié 


Les équations différentielles à argument dévié appartiennent à la classe des 
équations différentielles fonctionnelles (EDF). Elles décrivent l’évolution des 
variables dépendent des valeurs prises dans le passé ou le futur. On distingue 
trois catégories principales, les équations différentielles avancées, les équa- 


tions différentielles retardées et les équations différentielles mixtes. 


Equations différentielles avancées 


Sous leur forme la plus simple, ces équations s’écrivent comme suit : 


s (= f(t, z(t), e (t+ u)), a > 0 
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où f est donnée. Autrement dit, la valeur de la dérivée à un instant t de la 
variable liée x, ne dépend pas seulement de la valeur de x à l’instant t, mais 


aussi du futur. 


Equations différentielles à retard 


Sous leur forme la plus simple, ces EDR s’écrivent comme suit : 
a’ (t) = fa, x(t), s(t-r)), T> 0, 


où f est donnée. 
Autrement dit, la valeur de la dérivée à un instant ¢ de la variable liée x, ne 
dépend pas seulement de la valeur de x à l'instant t, mais aussi des valeurs 


prises avant l'instant t. 


Equations différentielles à argument dévié de type mixte (avancées- 


retardées) 


Sous leur forme la plus simple, ces équations s’écrivent comme suit : 
x (t)=f(t,c(t),c(t+p),x(t-7)), y, Tr >0, 


où f est donnée. Autrement dit, la valeur de la dérivée à un instant t de la 


variable liée x, dépend à la fois du passé et du futur. 


2.2.2 Equations différentielles à retard 


Les équations fonctionnelles à retards peuvent être classées comme linéaires 
ou non linéaires, autonomes ou non autonomes, périodiques ou non ou encore 
selon les types des retards. Dans ce paragraphe, on s'intéresse à donner une 
classification des EDR selon les types de retards cités dans la littérature où 
on distingue deux classes principales, la première est dite "d’équations diffé- 


rentielles retardées" et l’autre "d’ équations différentielles de type neutre". 
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Equations différentielles retardées 


Dans cette premiére classe, on trouve les sous catégories suivantes : 


` 


Equations différentielles à retard constant Sous leur forme la plus 
simple, les équations différentielles à retard constant, discret ou ponctuel 


s'écrivent comme suit : 
x (t) = f (t,x (t), x(t- 7)), 


où f est donnée et 7 est une constante positive. A titre d'exemple, on trouve 


ce type d’équations dans le modèle des mouches de Nicholson [A]. 


` 


Equations différentielles à retard variable Le retard dans ce cas va- 
riant dans la variable temporelle ou bien dépendant de l’état. A titre d’exemples, 
on trouve les équations différentielles 4 retard dépend du temps dans les mo- 
dèles de transport [49] et celle qui dépend de l’état dans le modèle décrivant 
l’évolution d’une population de poissons dont les larves consomment une 
nourriture [10] ou le problème de deux corps en électrodynamique classique 
[25] B&I [30]. 
Equations à retard variant dans le temps 


Sous leur forme la plus simple, ces équations s’écrivent comme suit : 
a (t)=f(t,c(@),2(¢-7(é))), 


où f est donnée. 
Equations à retard variable dépendant de l’état 


Sous leur forme la plus simple, ces équations s’écrivent comme suit : 


a(t) = f (t,e(t),2(t-7(#()))), 
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où f est donnée. 

Des fois cette classification n’est pas suffisante dans l’étude, ce qui conduit à 
ajouter des contraintes supplémentaires relatives au retard ou à sa dérivée et 
qui mène ensuite à définir de nouvelles sous catégories d'équations à retard 
variable comme : 

Equations à retard variable arbitraire : Le retard et sa dérivée ne sont 
pas limités. 

Equations à retard majoré : Cette sous-catégorie suppose la connaissance 


d’une valeur maximale sur le retard 
0 < T(t) < Tmax- 


Si T(t) = 7 est constant, il reste en pratique incertain et la contrainte ci- 
dessus assure un intervalle borné. Ce type de retard a été très largement 
abordé dans la littérature. 

Equations à retard (bi-borné) : Cette sous-catégorie est moins abordée 


que le cas précédent dont on suppose que le retard vérifie la contrainte 
Tmin < T (t) < T max- 


Equations à retard variant lentement dans le temps 


T (t) est une fonction dérivable presque partout telle que 
T(t) <A<1, 


qui indique alors une limitation sur la vitesse de variation du retard et que 
ce dernier varie lentement dans le temps autrement dit que les informations 
retardées arrivent dans l’ordre chronologique. 

Equations à retard variant modérément dans le temps 


T (t) est une fonction dérivable presque partout telle que 


T(t) < À avec À > 1. 
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Equations à retard variant rapidement dans le temps : Dans cette 


sous-catégorie, il n’y a aucune contrainte sur le retard et sa dérivée. 


Equations différentielles à retard distribué Ce type d’équations a été 
traité dans la littérature, pour la conception d’observateurs pour des systémes 
non linéaires ; pour prouver la stabilité des systèmes linéaires de type neutre ; 
la stabilisation robuste des systèmes neutres incertains et la commande ro- 


buste. Sous leur forme la plus simple, ces équations s’écrivent comme suit : 


v(t) = -as (t) -8 f 2(t-a)an(a, 


Par exemple, la chauve souris, à la chasse, étant aveugle, elle émet des sons, 
pour utiliser les parois des grottes, afin de localiser sa proie. L’écho obtenu 
par le rebondissement de ces cris représente le retard qui dépend de l’état, 
qui est le prédateur. On trouve ce type d’équations Aussi dans le modèle de 
sida ou le modèle de dynamique des populations présenté par Volterra 
en 1934 où il a utilisé un terme de retard distribué pour examiner un effet 


cumulatif sur le taux de mortalité d’une espèce. 


Equations différentielles à retard inconnu Dans ce cas, aucune hypo- 


thèse sur le retard n’est considérée qu’il soit constant, variable ou distribué. 


Equations différentielles à retard de type neutre 


Ces équations notées en abrégé EDN différencient des EDR. par le fait 
que la dérivée de l’état au temps actuel dépend non seulement des valeurs 
de l’état passé mais aussi de la dérivée d’ordre le plus élevé intervenant dans 


l’équation du temps passé. 
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Sous leur forme la plus simple, ces équations s’écrivent comme suit : 


d 
q Pe E-r 0) = S (te (t), a (t 7 (0), 


où f est donnée et D est un opérateur. Par exemple, on trouve ce type 


d’équations dans le modèle de réseaux distribués [49]. 


2.3 Résultats d’existence et d’unicité 


Dans cette section, on présente quelques résultats classiques sur l’exis- 


tence et l’unicité. 


2.3.1 Cas des équations fonctionnelles à retard 


Etant donné un nombre 7 > 0, C({a, b], R”), l’espace de Banach des fonctions 


continues définies sur [a, b] à valeurs dans R” muni de la norme de la conver- 


gence uniforme. Si [a,b] = [—7,0], on pose C = C(|[—7, 0], R”) et on désigne 


la norme d’un élément y € C par 


loll = sup ly]. 


P 
—T<t<0 


Si to ER, A > 0 et z € C ([to — T, to + A] , R”), alors pour t € [to, to + A], on 


définit x, E€ C par 
rls) = x(t +), 


pour tout s € [—7, 0]. 


Définition 2.1 [Ai] Si © est un sous ensemble de R x C, f : D — R” est 


une fonction donnée et "" représente ici la dérivée à droite, l’équation 


x(t) := f(t, x+), (2.1) 
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ou 


xls) := z(t + s), s € [—7,0], (2.2) 


est une équation différentielle fonctionnelle à retard sur Ð notée en abrégé 
EDR et le nombre 7 est appelé le retard. 

Il est clair que le cas 7 = 0 correspond au cas des équations différentielles 
ordinaires. Il est évident qu’une condition initiale appropriée au temps t = to 


exige la détermination de la fonction x sur tout l'intervalle [to — T, to], i.e. 


x(t) = Y(t), tE [to — 7, tol, (2.3) 


où Ÿ : [to — T, to] — R” est une fonction donnée supposée continue appelée 


la condition initiale de l’équation à retard (2.1). Ainsi, ’équation (2.1) peut 


se mettre sous la forme : 


r'():= fx), t> to 


(2.4) 
x(t) = y(t), t € [to — 7, tol, 


où 4 est une fonction donnée continue sur l'intervalle [to — 7, to]. 


Définition 2.2 Etant donné w € C et to € R, une solution de l'équation 


(2.1) est une fonction notée z(t) telle que x(t) = y(t) si t € [to — T; to] et 
satisfaisant (2.1) si t € [to, to + A] avec A > 0. Une telle fonction x(t) est dite 
solution de (2.1) à travers (4,4) et elle est notée souvent par 


x(t) = 2 (to,¥, f): (2.5) 
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Lemme 2.1 [4I] Etant données une fonction Y € C, to E R et f (t, Y) une 
fonction continue. La recherche d’une solution de l’équation à travers 


(to, Y) est équivalente à la résolution de l’équation intégrale 


Tto = Y, 
x(t) = ÿ(0) + Si f(u, £u )du, t > to. 


Théorème 2.1 (Existence) Pour l'équation (2.1). supposons que Q est 


un sous ensemble ouvert de R x C et f € C(Q,R") une application continue 


sur Q. Si (to, Y) € Q, alors il existe une solution de l'équation (2.1) passant 


par (to, y) i 


Définition 2.3 [41] On dit que la fonction f(t, y) est lipschitzienne par rap- 


port à y sur un compact K de R x C s’il existe une constante k > 0 telle 


que, pour tout (t,w;) € K, i=1,2,ona 


F(t, p1) — Ft, De) < k [br — dol. (2.6) 


Théorème 2.2 (Unicité) JJI} Supposons que Q est un sous ensemble ouvert 


de R xC, f : Q — R” est continue et f(t, w) est lipschitzienne par rapport 
aw sur tout sous ensemble compact de Q. Si (to, Y) € Q, alors il existe une 


solution unique de l’équation passant par (to, Y). 


2.3.2 Cas des équations fonctionnelles de type neutre 


Définition 2.4 [42] Supposons que 2 est un sous ensemble ouvert de R x C 


d’éléments (t, Y). Une fonction D : Q — R” [est dite atomique au point 3 


de Q si D est continue ainsi que sa première et seconde dérivée au sens de 
Fréchet par rapport à w et Dy, sa dérivée par rapport à w, est atomique en 


b de Q. 


1D west pas nécessairement linéaire. 
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Définition 2.5 [42] Supposons que 2 est un sous ensemble ouvert de R x C, 


D : Q — R” et f : Q — R” sont deux fonctions données continues avec D 


atomique en zéro. La relation 


d 
Bt Lt) = f(t, xt), (2.7) 


est dite une équation différentielle de type neutre notée en abrégé EDN. 


Théorème 2.3 (Existence) Si Q est un sous-ensemble ouvert de R x C 
et (to, Y) € Q, alors il existe une solution de l'équation passant par 


(to, Y). 


Théorème 2.4 [/2/ (Existence et unicité) Si Q est un sous-ensemble ouvert 


de R xC et f (t, Y) est lipschitzienne par rapport à ù sur tout sous-ensemble 


compact de Q, alors pour tout (to, Y) € Q, il existe une solution unique pour 


l'équation passant par (to, Y). 


2.4 Quelques méthodes de résolution des EDR 


2.4.1 Méthode des étapes 


La méthode des étapes (dite aussi "méthode pas à pas", "méthode des 
pas" ou "méthode des intégrations successives") permet de résoudre numé- 
riquement les EDR et les EDN et permet par la même occasion d’établir 
Vexistence et l’unicité de la solution. Elle fut présentée en 1965, par R. Bell- 
man pour des retards constants. D’autres comme El’sgol’ts et Norkin (1973) 
ont montré qu’elle restait aussi valable pour les retards variables, 4 condition 


que le retard ne s’annule jamais. 
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Pour fixer les idées, considérons l’équation fonctionnelle linéaire retardée avec 


coefficients variables suivante : 
x! (t) = a (t) x (t) + ag (t) x (t —7) pour tout t € [0,7] (2.8) 
ot) =o (4) pour tout t € [—7, 0]. 
Dans le cas où a; et a sont deux constantes réelles, alors l’équation est dite 
de type "Frisch-Holme" [36]. Maintenant, on va résoudre cette équation par 
la méthode des étapes. Le principe de cette méthode est de chercher des 
solutions sur des intervalles de type [kr, (k + 1) 7] où k € N en suivant les 
étapes suivantes : 
1° étape : Dans l’intervalle [-r, 0], la fonction x (t) est la fonction donnée 
y(t), donc l’équation est résolue dans l'intervalle [—7,0] et on désigne par 
xo (t) cette solution. Il faut signaler ici que si t € [0,7], alors t — T résidera 
dans [—7, 0]. 
2ème étape : Dans l'intervalle [0,7], si t € [0,7], alors t — 7 résidera dans 
[—7, 0], donc z (t — T) = xo (t — T) dans l'intervalle [0,7] et le système 


devient : 


8. 
se 7 
cH 
S” 

|l 


ay (t) x (t) + az (t) £o (t— 7T) pour tout t € [0,7] (2.9) 
qui est un problème à valeur initiale pour une équation différentielle ordinaire 
(EDO) où zo (t — T) = y(t — T) est connue. Ainsi, on résout cette EDO dans 
[0,7] en utilisant la condition initiale x (0) = y (0) et on désigne par x: (t) 
cette solution dans [0,7]. 


3%™° étape : Dans l'intervalle [r, 27], le système devient : 


w' (t) = a (t) x (t) + az (t) zı (t—7) pour tout t € [7,27] (2.10) 


Cette EDO avec la condition initiale x1 (T) est à son tour peut être résolue 


pour trouver la solution xə (t) € [7,27] et ainsi de suite. 
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Exemple 2.1 Considérons le cas particulier suivant : 


dg =ax(t—7T), pour t € [0,27] 
z(t) =p (t)=1 -Tr<t<0. 


étapes Dansi- 
xr(t)=1. 
2ème étape : Intégration dans [0,7] 


L’intégration des deux membres de L’EDR de 0 a t, donne 


v(j=af x(s — T)ds + x (0). 


Comme 0 < s < 7, alors -T < s—7 < 0. Sachant que x(t) = 1 pour 
t € [-7, 0|, alors 


x(s—7T)=l1, 


pour s € [0,7], ce qui conduit à 
z(t) = at +1, 


dans l'intervalle [0,7]. 
3°”¢ étape : Intégration dans [r, 27] 


L’intégration des deux membres de 7 à t, donne 


59 =a f a(s- 7)ds+2(7); 


Comme 7 < s < 27, alors 0 < s — T < T. Sachant que z(t) = at + 1 pour 
t € [0,7], alors 


x(s—T) = a(s—-7T) +1, 


pour s € [7,27], ce qui conduit à 
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dans l'intervalle [r, 27]. 


Finalement, on a obtenu 


t EDO avec condition initiale 
zw’ (t)=a 
0.7 (t) 
x(0)=1 
a (t)=o7%(t-—T) +a 
ag |) “=e 
x(T) = ar +1, 


et la solution est donnée par 


at +1 OSET 
x(t) = £2 a 
D + (@— a?r)t — a + or Sr TES r. 
Si on prend a = —1, 7 = 1 avec la même condition initiale ọ (t) = 1, on aura 
li = 0<Ł<1 
— — 2t + — 1<t<2 
2 2 


0 


Shen. 


ewww EERERERTS 


N'E=mmmmmni es 


Fic. 2.1 — Représentation graphique de l'équation x’ (t) 
p(t) = 1 (plein) et y (t) = —1 (pointillé) [69]. 


—x(t — 1) avec 
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Exemple 2.2 On considère le système suivant : 


x(t) = Ax(t) + Ba(t —7r)+ Max(t-7), p.p.t>0, 7 >0, 


avec la condition initiale 


(2.11) 


o € H! = H! ([-7,0],R"), -r<0<0. (2.12) 


1” cas : Si A = 0, B et M sont des matrices constantes de dimension n x n, 


la solution de 


r(t) =x 


est 


(0) + f,[Bz(s — r) + Mz(s —7)]ds, t > 0, (2.13) 


mais dans l'intervalle t € [0,7] 


Le changement de variable 0 = s — T et la condition (2.12) sur l'intégrale 


(2.13), nous donnent 


x(t) = 2(0) + [5 "[B6(0) + Mé(6)1d6. 


(2.14) 


Alors, l’utilisation de (2.13) de nouveau, on peut calculer x (t) dans Pinter- 


valle t € [7,27] et ainsi de suite. 


2ème cas : Si A # 0 dans (2.11), la solution est plus facile à calculer si on 


réécrit (2.11) sous la 


u 
dt 


[x(t) — Ma(t —7)] = Ax(t) + Ba(t—7), t> 0, 7 > 0, 


forme 


avec z, € C ([-7, 0], 


(2.15) 


R”) et la condition initiale xo(0) = ¢(0), € C ([-7, 0], 


R”), 


—T < 0 < 0. Alors il existe une solution unique x continue dans (—T, co) don- 


née par 


x(t) = Ma(t—r)+e“"[6(0)—C¢(—-7)]+ foe“) (AM+B)a(s—r)ds. (2.16) 
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On note aussi que la solution générale de (2.11)—(2.12) est pour t > 0 
_ „At t A(t—s) : 
a(t) = ep0) + Le [Bx(s —7) + Ma(s — 7) |ds. (2.17) 


Ceci peut être vérifié par dérivation directe de (2.16) et (2.17). 


2.4.2 Méthode de Runge Kutta 


Les méthodes de Runge Kutta font partie des méthodes d’analyse numérique 
à un pas utilisée pour la résolution des équations différentielles ordinaires. 


Une méthode de Runge Kutta de rang q est une méthode à un pas avec 


q 
P(t;, h, Ti) = So uk; 
j=l 
q 
kj = kiti) =f fe +c; h, £i + D) 1 < j < g{2.18) 
i=1 


OÙ Wj, Cj, Gj, sont des constantes. 


Elle est 
explicite si aji = 0, Vi >j 


semi-explicite si aj; = 0, Vi > j 
implicite dans les autres cas. 


Puis on a la relation suivante, entre le rang q et l’ordre p de la méthode 


Par exemple 
q= 1 on a la méthode d’ Euler, 
q = 2 on a la méthode d’Euler améliorée par exemple, 
q = 4 on a la méthode " classique " de Runge Kutta, 
q = 5 on a la méthode de Runge Kutta-Merson par exemple. 


La méthode que l’on a choisie est la méthode" classique" de Runge Kutta de 
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rang 4 et d’ordre 4 : 


Ti+1 — Ti } 2 (ka } 2ko } 2ks } ka), (2.19) 
avec 
ky = f (ti, xi) 
h 
2 2 
h h 
k3 = tit pat sh] 


Adaptation au probléme avec retard La méthode "classique" de Runge 
Kutta est extensible au cas des équations différentielles héréditaires |68]. En 


effet, considérons PEDR avec condition initiale suivante : 


x(t) = f (t, z(t), z(t- 7)) t <0, (2.21) 


x(t) = O(t) T<t<o. 


Prenons h tel que tT = mh, m € N*. On a alors (2.19) avec 


ky = de) 


h h 
ko = f ti + 5, ti + zki t 1 , 
2 2 para 
2 
h h 
k3 = f ti + 5, tit zko, T 1 ; 
2 2 i—m+— 
2 
ka = f(ti+ h, x; + hkz, mi mu). 


Les valeurs intermédiaires nécessaires sont calculées à l’aide d’une interpola- 


tion polynômiale de la solution entre les nœuds i — m et à — m +1. 


Remarque 2.1 On peut résoudre numériquement, les équations différen- 


tielles à retard en adaptant par exemple les solveurs : 
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DELSOL, RETARD et RADARS basés sur la méthode de Runge-Kutta 
(Fortran), dde23, ddesd et ddensd basés sur la méthode de Runge-Kutta 
(Matlab) et la fonction Lambert W basé sur la méthode des étapes (Maple). 


2.4.3 Méthode de Laplace 


On va résoudre l’équation de type Frisch-Holme suivante : 


y(t) =cy(t)+by(t—7), pourt>0. 
y(t) = h(t) Vt € [—7,0]. 


(2.22) 


Estimation exponentielle 


On va estimer la solution du probléme (2.22) par le théoréme suivant : 


Théorème 2.5 [2% Soit y(t) une solution de l'équation (2.29). Alors il 


existe deux constantes positives a et B telles que 


ly(t)| < a [|All exp (Bt), (2.23) 


où [h] = sup |h(0)]|. 
—T< 0<0 


Preuve. On intègre l'équation de type Frisch-Holme (2.22) de 0 à t et on 


[vus = e [uoas +6 f uis—r) as 


t—T 


y(t) = h(0) + ef y(s)ds + bf y(s) ds, 


=T 


arrive à 
donc 


ce qui implique que 


POESIE i lys) ds + [ol f ly(s)| ds, 
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i.e., 


all + (lel + 18) | lys) ds + [ol J mea 


e 
TRS 
~ 
>= 
A 


IA 


G + [ol 7) {All + (lel + hf ly(s)| ds. 


En utilisant le lemme de Gronwall, on obtient 


IA 


KI] (1+ Jol 7) [|All exp (Cel + 16/)t), 


< a|[hll exp (54), 


A 


avec a = 1 + |b|7, 8 = |c| + |b|. m 


L’équation caractéristique 


L’équation caractéristique d’une équation différentielle homogène linéaire à 
coefficients constants est obtenue à partir des solutions non triviales de la 


forme exp(At). L’équation 
y (t) = cy(t) + bylt = 7), 
possède une solution non triviale de la forme exp(At) si est seulement si 
Aexp(At) = cexp(At) + bexp( A(t — 7)). 


Ceci implique que 


À = c+ bexp(—7A), 


ou encore 


f(A) = À — c — bexp(—TÀ) =0. (2.24) 


f(A) = 0 est appelée l'équation caractéristique de (2.22). 
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La solution fondamentale 


Soit Y(t) une solution de l’équation (2.22) ayant pour condition initiale 


0 sit <0, 
h(t) = (2.25) 
1 sit=0. 


Y (t) est appelée solution fondamentale du problème (2.22). 


Théorème 2.6 La solution fondamentale Y (t) du problème est 


donnée par sa transformée de Laplace 


LY HA) = TO (2.26) 
Aussi pour y > B (8 = |c| + |b|) 
o f pAn O w T HT exp(At) 
Y(t) = iy Fay Oo m an FO A (2.27) 


Preuve. Puisque Y (t) satisfait l’estimation exponentielle (2.23) alors L(Y) 


existe (L(Y) est la transformée de Laplace de Y). On a 


i exp (—At) y'(t)dt = i exp (—At) c y(t)dt + [ exp (—At) b y(t — r)dt. 
a  — 


J 
(2.28) 


On intègre J par parties pour avoir 
J = [exp (—At) y(t)]6 + [> exp (—At) y(t)dt. 
On remplace par Y (t), on obtient 
J = [exp(-At) Y@]o + [> exp (—At) Y(t)dt, 


= -14 f dexp (=A) Y (t)dt. 
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L’équation (2.22) devient 
—1 +f Aexp(—At) Y(t)dt = | exp (— At) cY (t)dt + 
0 0 
+f exp (—At) cY (t — T)dt. 
0 


Donc, 


-14A\£L (Y(t) (A) = cL(Y (t))+ bf exp (—A(t + 7))Y (t)dt 


0 


= c£(Y (t)) (A) + bexp(—Ar) x J exp (—At) Y (t)dt 


—T 
oo 


+ exp (—At) vod 
) 


= cL(Y (t)) (A) + bexp(—Ar) L (Y (t)) (A) + 


+bexp (—Ar) j exp (—At) Y (t)dt, 


= cL(Y (t))(A)+ — (APD L(Y (t)) (À). 


(A — c — bexp (—Ar)) £(Y (t)) (A) = 1. 


On conclut que 


LY 0) A) = y= Co 7 Foy 


La preuve est achevée. m 


Représentation intégrale de la solution 


On va représenter la solution y(t) du problème (2.22) ayant pour condition 


initiale h(t) en fonction de la solution fondamentale Y (t). 


Théorème 2.7 La solution y(t) du problème peut étre représen- 


tée sous la forme 


y(t) =Y (t) h(0) + of y (t — 0 — 7) h(0)d6. 


—T 
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Preuve. On applique la transformation de Laplace à chaque terme du pro- 


blème (2.22), on obtient 


[ sud a E Gas sf OT ae 


a, amaai 
I 


(2.29) 


On intègre J par parties et on trouve 
I = [exp (-At) yel + | Xexp(—At) y(t)dt, 
0 


= —h(0)+ J reor uthat, 


0 


= —h(0) + AL (y) (A). 


L’équation (2.29) devient 
KOHLA] = è 


JAD exp (ate + 7)) vide 
y) (A) + bexp(—Ar) 

exp (—)t) y(#)dt + [ exp (—At) y(#)de| | 
y) (A) + bexp(àr)£ (4) (à) 
+bexp(—àr) i; exp (—At) y(#)dt. 


I 
Q 
X A à 
PS EL SN LES 
l © 


Donc, 
h (À) L(y) (A) = h (0) + bexp(-Ar) | exp (—At) h(t)dt, 


ceci implique que 


1 0 
LO) = ay h Gites [ex 107] | 
Donc, 
a ony) exp(—ArT m — 
wo = [TS [nor ES x] m M0) (By) dx, 


= y(t}h(0) - J Ds jeran x f exp (A) h(@)dd| a 
7 a 


Ti 
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Pour calculer J4. On définit une fonction w : [—T; +oo[— [0,1] telle que 
0 sifd>0, 
w (8) = 
1 si0 <0. 


= f EE 40 


T 


En posant s = T + 0, on obtient 


h = a exp(—As)h(s — T)w (s — T) ds, 
= L(h(=r+-)w (=r ++-)) (A). 


Maintenant, on va calculer 


PERAD ae 
f (A) 


f POAN s 
( 


REO) exp(At)L (Y) (A) £(A(-7 + -Jw (=r +-)) (A) A, 


exp(At)L£(Y x [h(—7 + -)w (—7 + -)]) (A) dà, 


Y (t — s) [h(—7 + s)w (—7 + s)| ds, 


Y (t — s) h(=7 + s)ds, 


= b y * [h(—7 + ju (—7 + -))) (©, 
i 
i 


car w (—T + s) = 0 dans [r,t]. Finalement (0 = s — T) 
y(t) = Y (t) A(0) + of y (t — 0 — T) h(0)d6. 


—T 


Ce qui achève la démonstration. m 


42 


Chapitre 2. Equations différentielles à retard 


2.5 Modèles d’équations différentielles à re- 


tard 


2.5.1 Equation de Nicholson (mouches-moutons) 


En 1950, le célèbre biologiste entomologiste australien Alexander J. Ni- 
cholson a mené une longue série d'expériences visant à en apprendre davan- 
tage sur des populations de mouches dévoreuses de viandes responsables de 
90% des myiased?| ovines qui menacent les élevages de plusieurs pays comme 
l'Australie, la Nouvelle Zélande et l Afrique du sud [65! 64 [66]. 

Cette mouche diptère ayant un corps rond à ovale de longueur varie de 
4.5 à 10 millimètres avec des yeux rougeâtres et un corps verdâtre ou vert 
bleuté avec des reflets cuivrés fait partie de la famille des "Calliphoridae" 
et elle est connue sous le nom "lucilie cuivrée australienne" ou tout sim- 
plement "mouches du mouton australien" ou en Latin "Lucilia cuprina" ou 
"Phaenicia cuprina". 

En outre, elle a deux paires d’ailes, la première paire étant des ailes mem- 
braneuses et la seconde paire étant des ailes postérieures réduites et modifiées 
connues sous le nom "d’haltéres" qui sont utilisées pour la stabilisation du 
vol. 

Le cycle de développement de cette mouche comprend quatre stades de crois- 
sance : ceuf , larve , pupe et adulte. La lucilie femelle gravide attirée par les 
plaies ou les replis laineux malodorants et humides des moutons pond en 
moyenne 250 œufs sur la peau de l’animal et qui vont éclore et se muent 
en larves carnivores aprés une période d’incubation ne dépasse pas 24h. Ces 


asticots se nourrissent des sécrétions des plaies et des tissus sous-jacents du 


?Une myiase désigne tout parasitisme d’un être vivant par les larves de diptéres para- 


sites. 
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Fic. 2.2 — La lucilie cuivrée. 


mouton pendant trois stades larvaires de durée de 4 à 5 jours. 

Après la phase larvaire, les larves complètement développées se laissent 
tomber et s’enfoncent dans le sol pour se transformer en pupes donnant des 
nouvelles mouches jeunes. 

Motivés par les données expérimentales obtenues par Nicholson, Gurney, 
Blythe et Nisbet ont proposé en 1980, l'équation retardée ci-dessous 
qui décrit l’évolution de la dynamique d’une population au cours du temps 


dN (t) 
dt 


= 6N (t — T) exp (=) -ôN (t), (2.30) 


où 
N (t) représente l'effectif de la population à l'instant t (les lucilies cuivrées 


adultes à l’instant t) 
dN (t) 


représente le taux d’évolution du nombre de la population 
B est le maximum de la croissance quotidienne d’ceufs par individu 


k est le nombre maximal d’individus que le milieu peut supporter] 


3La capacité de charge du milieu. 
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par les plaies ou les replis 
laineux malodorants et 
humides des moutons A 


Mouche femelle gravide attirée — Mouches adultes 


—s SN LG 
a 4.5-10mm > 
une odeur caractéristique 
souvent non perçue par 
l'homme ` Nouvelles 


y Mouches 
vy 
Environ lo oeufs pondus + 


TE 77 


vi Les pupes grandissent dans le 
Les larves complètement sol pour se transformer en 
développées se laissent nouvelle mouche 

tomber au sol 


Les larves se 
nourrissent de la chair 
du mouton 


Fic. 2.3 — Cycle de développement de la lucilie cuivrée. 


ô est le taux de mortalité par individu (jour~') 

T est la durée de la phase de maturation (le cycle de développement). 
Plusieurs généralisations et une recherche très active s’est amorcée récemment 
sur cette équation, par exemple en 2004, Liu et Ge [55], ont étudié l’équation 
différentielle de type neutre du premier ordre suivante : 


“ [x (t) — ex (t — T (t))| = —a (t) x (t) + b(t) æ (t — r (t)) CPOE), 


En conservant la même dénomination; "équation de Nicholson", de nou- 
veaux modifications du modèle ont été largement considérées dans l’étude 


des croissances d’autres espèces. 


2.5.2 Problème de deux corps en électrodynamique 


classique 


Dans (voir aussi [43]), Driver a étudié un problème de deux 


corps en électrodynamique classique avec une interaction retardée. Il a consi- 
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déré le mouvement de deux particules chargées se déplaçant le long de l’axe 
des x et substitué l’expression pour le champ d’une charge mobile, calculée 
à partir du potentiel de Liénard- Wiechert, dans la loi de la force de Lorentz- 
Abraham. La réaction de rayonnement est omise, mais le retard est inclus en 
raison de la vitesse limitée de la propagation, c, des effets électriques. la mo- 
délisation mathématique de ce problème conduit à un système d’équations 
différentielles à retard. En partant de ce modèle, il a obtenu un système de 
six équations différentielles décrivant l’évolution des états, des vitesses et des 
retards temporels. 

Pour décrire son modèle, nous désignons par x; (t) (i = 1,2) les positions 
des deux charges ponctuelles sur l’axe dans un système inertiel donné à l’ins- 
tant t. Soit vu; (t) = x; (t) (i = 1,2), les vitesses des charges. Comme déjà 
mentionné, nous omettons la réaction de rayonnement mais nous permettons 
un champ électrique externe, Eert (t, £), dans la direction x, qui est supposé 
continu sur un ensemble ouvert D; dans le plan (t,x). Alors, l’équation du 
mouvement d’une particule chargée 7 est donnée par : 

mv; (t) 


a-ge A 


Ej (t, xi (t)) + qiEext (t, xi (t)), i,j € {1,2}, j # i, 
(2.31) 
où m; est la masse au repos et q; est la grandeur de la charge i, c est la 
vitesse de la lumière dans le vide, et E; (t, x) est le champ électrique à (t, x) 
dû à une autre charge 7 # i. Le champ magnétique de la charge j n’est pas 
impliqué dans ce cas unidimensionnel. 
Le champ à l'instant t et au point x produit par la charge j est supposé 
calculé à partir des potentiels de Liénard- Wiechert. L’expression de ce champ 
implique un décalage temporel t — Tji, représentant l'instant où un signal 


lumineux devrait quitter la charge j pour arriver à x; (t) à l'instant t. Donc, 
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le retard doit être une solution de l’équation fonctionnelle 
Tji (t) = |x: Œ) — z; (t — Ty (t))| 7e. (2.32) 


Il est clair que 7;; (t) ne peut pas être écrit explicitement. 

L'introduction des retards dans l'équation (2.31), nécessite de déterminer 
des trajectoires initiales des deux charges sur un certain intervalle approprié 
la, to]. 

Considérons maintenant ces trajectoires initiales et leurs extensions (2; (t) , x2 (t)) 
définies sur un certain intervalle [a, 3), où 8 > to, tel que 

(a) chaque x! (t) est continue et |x; (t)| < c pour tout t € [a, 8); 

(b) x2 (t) > x2 (t) et (t,x; (t)) € D pour tout t € to, 8); 

(c) les deux équations fonctionnelles 7, = |i (to) — z; (to — To) | /c ad- 
met des solutions 79,, i A j, i, j € {1,2}. 

Puis Driver prouve que (xı (t) , z2 (t)) est une solution de (2.31)— (2.32) 
si et seulement si elle satisfait au système de six équations différentielles à 


retard suivant pour t € (to, 3) : 


2 (8) = où (8), | 
(Div: (t) = (Dwy (t= TH (0) 


Ti QE c— (—1} v; (t — Tj (t)) | 
v! (t) 2 (—1)’ aic c+ (—1)' 9; (t — Tj (t)) + qi Eert (t, zi (t)) /mi 


(L—v2(t) /ep? TE) ce (1u (t ryj lt) 


(3 


(2.33) 

Où 754 (to) = Thi, ai = q1q2/ (Am €o Mi È) est une constante, et en particulier, 
€ est la constante diélectrique de l’espace libre et (i, j) = (1,2) ou (2,1). 

Driver [28] a montré que si les trajectoires initiales données satisfont à 

la condition (a) pour a < t < to, la condition (b) en to, et la condition (c), 

et si Eert(t, x) est lipschitzienne par rapport à x dans chaque sous-ensemble 


compact de D; et si la vitesse initiale de chaque particule est lipschitzienne, 
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alors il existe une solution unique. Cette solution peut être continuée d’exister 
tant que les charges n’entrent pas en collision (lim z; (t) = lim x (t) quand t 
se rapproche de l’extrémité droite de l’intervalle maximal d’existence) et ni 
(t, xı (t)) ni (t, £2 (t)) ne s'approche de la frontière de D4. Il faut signaler ici 
que dans l’article de Driver et Norris [BI], la condition de Lipschitz ci-dessus 
pour les vitesses initiales est relâchée à l’intégrabilité de la vitesse initiale sur 
[a, to]. 

Dans Driver [24], un cas particulier a été donné où les positions et les vi- 
tesses des particules à un instant donné détermineront l’état du système. 
Plus précisément, dans cet exemple d'équations électrodynamiques du mou- 
vement, les valeurs instantanées des positions et des vitesses des particules 
détermineront leurs trajectoires, si les solutions sont définies pour tous les 
temps postérieurs. Cette propriété était fréquemment conjecturée, affirmée, 
ou implicitement supposée, comme dans la mécanique newtonienne et comme 
indiqué par la longue liste de références citeés dans l’article de Driver [24], 


mais cette propriété ne doit pas être attendu pour les équations électrody- 


namiques générales. Dans le cas où Eer (t,£) = 0 pour tout (t,x) € R? et 


4192 > 0 (les deux charges sont de même signe), lim [to (t) — x, (t)| = co et 
lvi (t)| < c pour tout t > a. C’est un résultat très intéressant qui prouve que 
le retard Tj; (t) peut devenir non borné, tel qu’on obtient un système d’équa- 
tions différentielles fonctionnelles avec des retards non bornés dépendants de 
l’état. 

Dans le cas où les mouvements sont considérés tridimensionnels, alors on 
obtient un système d’équations différentielles fonctionnelles de type neutre 
avec retards dépendants de l’états, et le taux de changement de v; au temps 


actuel dépend aussi de sa valeur historique vi (t — T;;). Plus précisément, si 
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nous introduisons un vecteur unitaire 


Li — Tj (t — TH 
y= HO TH) 
CT ji 


et une quantité scalaire 


1 
Vij = 1 = 78 (t — Tyi) ui, 


comme dans [27], le point ”.” désigne le produit scalaire dans R? (sans oublier 


que dans ce cas, £1, £2 sont des vecteurs dans R?), alors la loi de force de 


Lorentz donne 


/ 
w(t) — 2 Ce LE 


[Ej + (ui/c.E;) (u; — vi/e) — (vife-us) E;], 
(2.34) 


Mi 


où E; est le champ électrique retardé (valeur vectorielle) arrivant à x; à 


l'instant t de la particule j. Ce champ, dans RÌ, peut être trouvé à partir des 


potentiels de Liénard- Weichert tels que 


keqi 
B= a a [us — oi (¢— e =a] 
ji lij 
kqi 
bi mxiet (2-35) 
Tji Vij 


où k > 0 est une constante en fonction des unités, et x désigne le produit 


vectoriel croisé dans R°. L’adaptation dynamique de 7;; est donnée par 


26-22 ae (2.36) 


ji 
CY 45 


D’aprés ce qui précéde, le mouvement de chaque particule est influencé par 
les champs électromagnétiques des autres, et en raison de la vitesse finie de la 
propagation de ces champs, les équations du modèle décrivant le mouvement 
des particules chargées via une action à une distance impliquera des retards 


temporels qui dépendent de l’état du système entier. 
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Driver et Hoag [26} [44] ont noté que si l’on considère que les lois de base 
de la physique sont symétriques par rapport au renversement du temps, donc 
l’existence de ces retards implique l'existence des termes avancées dans les 
équations, et donc on est conduit à un système d’équations différentielles 
fonctionnelles avec des arguments mixtes [44] ou de type neutre [26]. 

L’héritage des travaux de Driver et ses collaborateurs en électrodyna- 
mique est extraordinaire. Malgré qu’ils ont été publiés depuis 55 ans, ils 
marquent la naissance de la théorie des équations différentielles à retard dé- 
pendant de l’état et ils ont été couronnés de nombreux succès dans l’inspira- 


tion pour le développement théorique et pratique pour de nouveaux résultats. 


Cependant en pratique, beaucoup d’EDR peuvent déduites à partir d’équa- 
tions aux dérivées partielles (EDP) ou d’équations aux dérivées partielles à 
retard (EDPR). On illustre cela avec trois exemples : l’évolution au cours 
du temps des populations de cellules souches hématopoïétiques, une ligne de 


transmission sans pertes et un modèle de réaction-diffusion. 


2.5.3 Hématopoïèse 


Il ne s’agit pas ici de faire une présentation biologique détaillée de l’hé- 
matopoïèse, mais plutôt de donner une idée de ce qu’est l’hématopoïèse, 
de présenter un modèle mathématique pionnier et de citer quelques travaux 
dédiés à cette thématique et qui nous paraissent incontournables. L’hémato- 
poïèse est un mot venant du grecque haimatôpoiein, d’après poiein "faire" 
et plus précisément, c’est la production et le renouvellement quotidien d’une 
quantité très importante de cellules sanguines (les globules rouges lors de 
l’érythropoïèse, les globules blancs lors de la leucopoïèse et les plaquettes 


lors de la thrombocytopoïèse) dans les organes hématopoïétiques (moelle 
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osseuse chez l’adulte, foie et rate chez l’embryon). 

Cette production des éléments figurés du sang est assurée par les cel- 
lules souches hématopoïétiques (CSH) dites primitives et qui ont la capacité 
d’auto-renouvellement et de différenciation vers l’une ou l’autre des lignées 
cellulaires. Elles se différencient en quatre sous-populations selon leurs fonc- 
tions et leur état de maturation : cellules souches (immatures), progéniteurs, 
précurseurs (cellules en cours de maturation.) et cellules matures (différen- 


ciées). 


== yi 
e »,Auto-renouv ellement Moelle 


Sang 


x2 Osseuse 


Cellules 


souches 
we i e 


Contrôle à court terme 


1 
i 
i B 


Cycle cellulaire Progéniteurs 


Taux de prolifération 
Apoptose 


Durée de prolifération 
Apoptose 


Auto-renouvellement 

Différenciation Maturation A, B : Compartiments de prolifération 

Apoptose C, D : Compartiments des cellules 
différenciées 


Apoptose 


Fic. 2.4 — Schéma représentant les différentes étapes de l’hématopoïèse 


Les premiers modèles mathématiques remontent à la fin des années 70 et 
parmi les travaux pionniers qui forment les bases de l’étude mathématique 
de l’hématopoïèse, on trouve les études de M.C. Mackey en 1978 [59] sur 


l’hématopoïèse périodique et l’anémie aplasiqud!| 


4Anémie due à un dysfonctionnement du tissu médullaire, incapable de produire cor- 


rectement les éléments figurés du sang. 
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Nous allons maintenant présenter le premier modèle de Mackey [59] qui 
est inspiré par des travaux antérieurs de Lajtha en 1959, Burns et Tannock 
en 1970 et qui décrit l’évolution au cours du temps de deux populations de 
cellules souches hématopoïétiques immatures, une proliférante (les cellules 
sont en phases G1, S, G2, M) et l’autre quiescente (les cellules sont en phase 
de repos Go). En partant d’un système couplé d’EDP hyperboliques du pre- 
mier ordre et grâce à la méthode des caractéristiques, il trouva un système 


de deux équations différentielles à retard. 


Fic. 2.5 — Le cycle cellulaire [2]. 


On a les équations de transport suivantes : 


? n(t,a) + En a) = (5 + B(N(D)n (a) a>0,t>0 


Y 


Z plt a) + <p(t,a) = -p (t0) O0O<a<T,t>0 
n (t,0) = 2p (t, T) t>0 
p(t,0) = J B (N (t))n(t,a) da t >0, 

(2.37) 
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avec les conditions initiales : 


n(0,a) = no (a) € L! ((0, +00) ,R*) 


(2.38) 
p (0, a) = po (a) € £1((0,7) ,R*), 
et la condition naturelle : 
im 7 a) = 0, 220, 
où 
a: l’âge qui est le temps passé par une cellule dans une des 
deux phases, quiescence ou prolifération. 
n(t,a): la densité des CSH en phase de repos à l’instant 
t > 0 et d'âge a > 0. 
p(t,a): la densité des CSH en phase de prolifération a 
Vinstant t > 0 et d'âge 0 <a<r. 
N (t): la densité de la population totale des cellules dans la phase de repos 
où N (t) = ae) (t, a) da. 
ô : Le taux de mortalité des cellules au repos (jour™t). 
y: le taux d’apoptose (mort naturelle, suicide cellulaire, ou mort 


cellulaire programmée) (jour” !). 
b: le taux de réintroduction (jour™t). 
la durée moyenne de la phase de prolifération, (la durée du 
cycle cellulaire) (jour). 
Dans le système (2.37), la premiére condition frontiére signifie que toutes les 
cellules proliférantes entrent immédiatement aprés leur division dans la phase 
Go. La deuxiéme condition frontiére décrit le flux des cellules quiescentes qui 
entrent dans la phase de prolifération. 


En vertu de la méthode des caractéristiques, on arrive à 


e Yp(0,Tr—-t) O<t<r 
Ce = (2.39) 
eVp(t—7,0) t>0. 
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En posant 
P (t) = f p(t,a) da, (2.40) 
0 
la densité des cellules dans la phase de prolifération à l’instant t > 0 et en 
tenant compte des conditions frontières et de (2.39), l’intégration du sys- 


tème (2.37) par rapport à l’âge donne le système couplé de deux équations 


différentielles à retard temporel discret suivant : 


N'(t)=—-(6+ B8(N(t))) N(t)+2e 76 (N(t—7))N(t-—7) t>7 


P@=—-yPO)+PNEN@)-e VPN ft) Nir) 427, 
(2.41) 


avec des conditions initiales pour t € [0,7]. 

Dans la première équation, la perte de la non prolifération des cellules à la 
phase de prolifération et à la différentiation est représentée par le premier 
terme de droite tandis que le deuxième terme représente la phase de produc- 
tion des cellules dans la phase de repos G à partir de la prolifération des 
cellules souches. Le facteur 2 compte pour l’amplification des effets de la di- 
vision cellulaire tandis que e~?" compte pour l’atténuation dans la phase de 
prolifération en raison de l’apoptose. Ce modèle a été modifié de nombreuses 


fois. 


2.5.4 Ligne de transmission sans pertes 


Maintenant, en vertu de la transformation de d’Alembert, on va voir com- 
ment peut on réécrire une EDP hyperbolique en tant que système à retard 
de type neutre [70] similaire à [49]. 

Cet exemple traite une ligne de transmission sans pertes, à une extrémité 


(x = 0), il y a une source de tension constante F et à l’autre extrémité 
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(x = 1), une capacité est branchée en parallèle avec une diode tunnef| 


R x=0 x=l 


FIG. 2.6 — Circuit de la ligne de transmission [70] 


Le courant i(.,.) et la tension v(.,.) sont des fonctions de t et x qui satisfont 
Véquation du télégraphe, donc on obtient les deux équations aux dérivées 


partielles du premier ordre de type hyperbolique suivantes : 


ði Ov Ov ði 
se SOEs pana — = <a< ; 
La az Dete a +a. eet Uae = 4, (2.42) 


avec les conditions aux limites 


E = Ri(t,0) + v(t,0) et i(t,l) = C, ot 1) + ig(v(t,2)), (2.43) 


où l est la longueur de la ligne, L et C sont l’inductance et la capacité du 
conducteur par unité de longueur, À est la résistance à l’entrée, C est la 
capacité en parallèle avec la diode tunnel, ig(v) est la courbe caractéristique 
courant-tension de la diode (une fonction polynômiale non linéaire) donnée. 


Les points d'équilibre (vo, to) satisfont une équation linéaire donnée par 


E = vo + Rio, to = ld (vo). 


5Ce genre de diodes est utilisé dans les amplificateurs haute fréquence de circuits os- 


cillants électroniques par exemple. 


99 


Chapitre 2. Equations différentielles à retard 


-0.0004 


Fic. 2.7 — Courbe caractéristique courant-tension et ligne de charge de la 


diode tunnel. [70] 


On suppose que le point de fonctionnement est donné par l'approximation 
du premier ordre 

i(va) ~ i(vo) +m (va = vo) ; 
où m est une constante positive. La formule de D’Alembert et les conditions 
aux limites (2.43) nous permettent d’obtenir l’équation différentielle linéaire 


de type neutre suivante : 


TEE) = Cage E-7) = AVE + BU r) (2.44) 


nya p Weee) ve 
a (AE) C, j i 


avec A = 


T = 2V LCI et E =: t — bl. 
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Ce qui prouve une transformation] du systéme de deux EDP de type hyper- 
bolique (2.43) en une équation a retard de type neutre. 


2.5.5 Equation de réaction-diffusion avec argument dé- 
vié 

L. Sergeeva [76], A.M. Samoilink et collab ont établi les conditions 
dans lesquelles il est possible de construire des solutions globales d’une EDP 
avec argument dévié dans la variable temporel. Ils ont montré l’existence et 
Vunicité de ces solutions en vertu du théoréme de point fixe de Banach et 
ils ont donné la structure de ces solutions et l’algorithme utilisée pour leur 
construction par l’utilisation de la méthode de séparation des variables et la 
méthode des suites récurrentes. 

Considérons le probléme aux dérivées partielles linéaire non homogéne 


avec argument dévié suivant : 


us (x,t) = p (t) Use (x,t + u) +q (x,t), (x,t) €Q, 
u (0,t) = u (l,t) = 0, te 


(2.45) 


A 


? 


avec des conditions aux limites de Dirichlet homogènes 


u (0,t)=u(l,t)=0, TER, 


où Q = {(x,t):0<x <l, te R} et u € R\ {0}. 


Remarque 2.2 Dans le cas où t > 0 (la solution n’existe en général que 
dans le futur), PEDP avec argument dévié (2.45) est de type parabolique et 


"4 


dite "équation de réaction-diffusion avec argument dévié" ou "équation de 


la chaleur non homogène avec argument dévié" [où u (x,t) est une fonction 


6La transformation n’est pas unique. 
TLe premier exemple a été étudié par Fourier en 1822 dans son œuvre révolutionnaire 


"Théorie analytique de la chaleur". 
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définie pour la position x et le temps t et représente les individus qui peuvent 
être différents d’un problème à l’autre ; ils peuvent représenter des concen- 
trations, des densités, des charges électriques, des températures, etc. La pre- 
mière dérivée par rapport au temps traduit la vitesse de formation, q (x,t) 
est une fonction appelée "la réaction", "le terme réactif" ou "la source" et 
elle découle de toute interaction entre les individus, p (t) est le coefficient de 
diffusion et p (t) Uze (x,t + u) est le terme de diffusion qui peut représenter la 
diffusion moléculaire ou un certain mouvement aléatoire d'individus dans une 
population et u est l’argument dévié. Les conditions aux limites homogènes 


de Dirichlet signifient qu’il n’y a pas d'individus sur la frontière. 


Pour simplifier l'EDP avec argument dévié (2.45), Samoilenko et all [73 
[72] [76], ont utilisé la méthode de séparation des variables sur le problème 
homogène qui a comme fonctions propres, 

krz 
Xx (x) = sin z 
associées aux valeurs propres 
rk\° 
Àk = (=) ,k =1,...,n pour certain n > 1. 
Il ont construit une solution globale u (x,t) du problème (2.45) sous la forme 


u(x,t) = yT (t) sin ave, (x,t) EQ. (2.46) 


L. Sergeeva [76] a supposé aussi que le terme source q (x,t) peut être re- 


présenté comme une somme des n premiers termes de la série de Fourier 


impaire 
= k 
q (x,t) = Sax (t) sin a (2.47) 
k=l 
ou 
1 i 
2 _ KT 
ak (t) = fa (£, t) sin Ty dE: (2.48) 


0 
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Donc, la substitution de (2.47) et (2.48) dans (2.45) donne l’équation à ar- 


gument dévié suivante : 


T! (t) + Aup (t) Tk (t + u) — qu (t) = 0, k=1.n (2.49) 


2.6 Récapitulatif 


Revenons un instant sur les principaux points abordés dans ce chapitre. 
Premièrement, on a présenté un bref historique sur les équations fonction- 
nelles à retard suivi par une classification des différents types de ces équations 
cités dans la littérature. Ensuite, on a proposé un rapide survol sur la théo- 
rie d’existence et d’unicité des EDR et des EDN et quelques méthodes de 
résolution comme la méthode des étapes, celle de Runge-Kuta et la transfor- 
mation de Laplace. Finalement, on a exposé cing modèles d’équations avec 
retard, les deux premiers modèles (le modèle de Nicholson et le modèle de 
deux corps en électrodynamique classique) peuvent être modélisés mathé- 
matiquement par des équations différentielles à retard discret ou variable ou 
bien par équations différentielles de type neutre tandis que le troisième et 
le quatrième phénomènes sont modélisés par des EDP sans retard mais leur 
étude mathématique conduit à obtenir des EDR ou des EDN. En effet, dans 
le modèle décrivant l’évolution au cours du temps des populations de cellules 
souches hématopoïétiques, la méthode des caractéristiques conduit à obtenir 
un système couplé de deux équations différentielles à retard temporel discret 
et dans le modèle d’une ligne de transmission sans pertes, la transformation 
de d’Alembert nous permet de réécrire l’'EDP hyperbolique en tant qu’une 
équation différentielle à retard de type neutre. 

Enfin, on a présenté un modèle d’une équation aux dérivées partielles 


linéaire non homogène à argument dévié où la méthode de séparation des 
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variables sur le problème homogène conduit à une équation différentielle à 
argument dévié. 

Pour plus de détails nous référons aux livres de Hale [41], Hale et Lunel 
[42] et le lecteur pourra consulter aussi le livre de Kolmanovskii et Nosov 


pour des exemples de modèles comportant des termes à retard. 
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CHAPITRE 3 


| 


Existence et positivité des solutions périodiques d’une 


classe d’équations différentielles fonctionnelles de type 


neutre du second ordre 


Sommaire 
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C chapitre est consacré à prouver l'existence de solutions périodiques 
positives d’une classe d’équations différentielles du second ordre non 
linéaires de type neutre. L’étude est faite en vertu d’une approche basée sur 
Vapplication du théorème de Krasnoselskii combiné avec quelques propriétés 


de la fonction de Green d’une EDO linéaire du second ordre. 
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3.1 Introduction 


Dans ce chapitre, nous sommes essentiellement intéressés par l’étude de 
l'existence de solutions périodiques positives d’une classe d'équations diffé- 
rentielles non linéaires du second ordre de type neutre qui apparaissent dans 
la modélisation de nombreux phénomènes rencontrés en physique, biologie, 
mécanique, économie, etc. 

On peut considérer que l’étude des équations fonctionnelles de type neutre 
commence avec les trois articles] (presque) oubliés de I.P. Kabakov (1946), 
Sokov (1946), I.P. Kabakov et A.A. Sokov (1946) sur les phénomènes de pro- 
pagation. À partir des années 60, ces équations ont été l’objet de plusieurs 
travaux. Notamment, on trouve les travaux de A. Bolinia et Myshkis (1960), 
Brayton et Mille -Luxriante (1964), Cooke et Krumme (1968), Cooke (1970), 
Halanay et Rasvan (1970) et Filimonor (1981), etc. 

Durant ces deux dernières décennies, de nombreux travaux ont établi plu- 
sieurs résultats sur l’étude des solutions périodiques. Certains mathémati- 
ciens ont utilisé des transformations pour se ramener à une équation plus 
simple ou à un système d’équations. D’autres chercheurs ont donné la solu- 
tion sous forme de série convergente vers la solution exacte et certains d’entre 
eux ont traité les équations fonctionnelles du second ordre de type neutre par 
des méthodes numériques telles que la méthode de Ritz, la méthode des diffé- 
rences finies, la méthode des éléments finis, la méthode des splines cubiques 
et la méthode des multi-dérivées. 

Dans ce chapitre, ces méthodes habituelles peuvent sembler inefficaces pour 


traiter l’équation différentielle non linéaire du second ordre de type neutre 


'Les trois articles ont été écrits probablement plusieurs années plus tôt mais publiés en 


1946 . 
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suivante : 
Sui) +P aA + ar) = Stel- rt) (3.1) 


+f (t, hı (x(t), (w(t — T(t)))), 


où p et q sont des fonctions positives continues. La fonction Q : R x R — R 


est différentiable, h1, gı : R — R et f : R x R x R — R sont des fonctions 


continues par rapport à leurs arguments. 

On va utiliser le théorème de point fixe de Krasnoselskii combiné avec quelques 
propriétés d’une fonctions de Green pour établir les résultats désirés. L’idée 
de l’approche est de convertir l'équation considérée en une équation intégrale 
qui s’adapte à l'application du théorème de point fixe de Krasnoselskii pour 


prouver l'existence de solutions périodiques positives. 


3.2 Préliminaires 


Dans ce paragraphe on donne quelques définitions et certains résultats essen- 
tiels qui représentent un outil important pour le reste de ce chapitre. 


Soit 


Pr = {a (t) €C(R,R): x(t+T)=x(t), T >0}, 


muni de la norme 


lællo = sup |a(t)| = sup |a(t)]. 
teR te [0,7] 
Alors (Pr, ||.||.o) est un espace de Banach. 


Puisque nous sommes intéressés par l’existence des solutions périodiques de 


l'équation (3.1), il est naturel de supposer que 


plt +T) = po), QU+T) =a etr(t+T)= 76), (32) 
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où 7 est une fonction continue, et T(t) > 7* > 0. En outre, nous supposons 
que 
T T 
| p(s)ds > 0, 1 q(s)ds > 0. (3.3) 
0 0 


Nous supposons aussi que les fonctions Q(t, x) et f(t, x, y) sont T—périodiques 


en t, ce qui est équivalent a 
Q(t +T,2) = Q(t,x) et f(t+T,2,y) = f(t,2,y). (3.4) 
D’abord, on considère l’équation suivante : 
x'(t)+a(t)x(t) = AM(x(t))+mt(t), (3.5) 
où a est une fonction continue 7 —périodique et non négative avec 


T 
| a(s)ds > 0, 
0 


m est une fonction continue T'—périodique, À > 0 est un paramètre et 


M : Rt — R* est une fonction continue. 


Lemme 3.1 Supposons que 


M 
lim 2 
xz—+o0 T 


= N € (0, +00]. (3.6) 


Alors, admet au moins une solution T—périodique positive si €, < E 
et A € (Ei, E2], où &, et £ sont définies comme suit : 


7 2 [exp (Jor a(u)du) — 1 be Ro [exp (le alu)du) — 1 
a NoT LL QT | 


avec 


Ro = [RE o= oxp (=f ttd), 


Qı = max V (x), V (x) = 
ze[0,Ro(o+1/0)] M(0) six € (—oo, 0) 


R(t) = [va (t,s)m(s)ds>0, Gi(t,s)= 
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Lemme 3.2 Supposons que et sont vérifiées et 
Ri [exp ta p(u)du) — 1 


QT 21, (3.7) 
où : 
tr exp (| plu) 
1 t€[0,T] f a (f noa) E M ; 
0 
et 


2 


Qi = (1 + exp (f veau) RÈ. 


Alors, il existe des fonctions continues T— périodiques a et b telles que 


b(t) > 0, fp a(u)du > 0 et 


a(t) + b(t) = p(t), —b(t) + a(t)b(t) = q(t) pour toutt € R. (3.8) 


Preuve. De (3.8), on obtient 
b' (t) + p(t)b(t) = 0? (t) + a(t). (3.9) 


Soit M(x) = x?, m(t) = q(t) et À = 1. Alors 


Par le Lemme si 


0<A=1< OT | 
alors, admet au moins une solution positive T —périodique continue b(t). 
Maintenant, comme 7 
| q(s)ds > 0, 
0 
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et 


il vient 


or = | La > 0. 


D'où le résultat désiré. m 


Lemme 3.3 [86] Supposons que les conditions du Lemme [3.1 sont vérifiées 


et ® € Pr, alors l’équation 


er | p(t) Gt) + g(t)x(t) = S(t), (3.10) 


admet une solution T—périodique. En outre, toute solution de s'écrit 


sous la forme : 
+T 
a= f G(t, s)®(s)ds, (3.11) 
JP exp |f; b(v)du + f? a(v)dv] du + 1 exp Bhs b(v)dv + E” a(v)do| du 


[exp (0 a(u)du) — 1 [exp (hs b(u)du) — 1 
(3.12) 


G(t,s) = 


Preuve. Définissons E, et Ep comme suit : 


Ey = exp (f aoan) — 1, E= exp ([ rear) =. 


Par un calcul direct, nous pouvons voir que (3.11) est une solution T —périodique 


de (3.10). 


Maintenant, on suppose que x(t) est une solution T—périodique de (3.10), 
en vertu du Lemme nous avons 


d? d 
PAU) +a(t)—=zx(t) + 


d d 
z FOLO +E) + a(D UD (E) = 20), 


dt 
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ce qui conduit à 


= LE n ne. f o ya) f “ew (/ Í o(u)du) ojiva: 
b(u)d Jr ( f | o(u)du) B(v)ds 
EE af" an “he a (fu jau) sp ( / “a(u)du) B(v)ds 
=- 55 zz © "a aa odo | l [ wae) du 


ey B(s ve z q eef aap) ra 
/ CE 8) 8(s)ds. 


Ce qui établit le résultat désiré. m 


Corollaire 3.1 La fonction de Green G définie par satisfait les 
propriétés suivantes : 

G@,t+T)=G,t), GU+HT,s+T) = G(t, s), 

: ae f saa 


—G (t, s) = a(s)G(t, s) — 
ðs 2 
exp( | b(v) dv) — 1 


: E exp( f ‘wee 
exp( / a (v) dv) — 


? 
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Remarque 3.1 Par un calcul simple, on obtient que la dérivée seconde par 


rapport à s de la fonction de Green G(t, s) est comme suit : 
exp( | b (v) dv) 
F : 
exp( | b(v) du) — 1 
0 


Le Lemme suivant est important pour établir l’existence. 


Lemme 3.4 Supposons que les conditions (3.2-(5.4) et sont satis- 


faites. Si x € Pr, alors x est une solution de l'équation si et seulement 


a(t) = Qtat— 7) — | p(s) E (t, 5) Q (s, a((s — T (s))) ds 


+] Git, s) [f (s, hı (x (s)), gı (x (s — 7 (s)))) (3.13) 
+ (a(s) + a°(s)) Q(s,2((s —T (s)))| ds, 


ou 


E(t,s) = exp( f; b (v) dv) 


7 exp([> b (v) dv) —1 (814 


Preuve. 


Ligne directrice : Par l’utilisation du Lemme 3.3 et certaines propriétés 
de la fonction de Green citées dans le Corollaire 3.1, on peut convertir 
l'équation (3.1) en l'équation intégrale (3.13). 
Soit x € Pr une solution de l'équation (3.1). D’après le Lemme on 
a 


rt) = / G(t, s) [20 (s, x((s — 7 (s))) (3.15) 


+f (s, ha (x(s)), gı (x (s — 7 (s))))] ds. 
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Une intégration par parties nous donne 


AOC 0) 


i G(t, 370 (s,x((s — T(s))) ds = -f © alt 370 (s,x((s — T(s))) ds. 


Posons 


Par le Corollaire [3.1] 


exp( | oidi ue 


I=- |Q(s,a((s—7(s))) | a(s)G(t, s) — 


t 


69 


Chapitre 3. Existence et positivité des solutions périodiques d'une classe 
d'équations différentielles fonctionnelles de type neutre du second ordre 


donc 
lS. GTa T =r) 
exp( | b(v) dv) 
alt +T)G(t,t +T) = 
x | a(t +T)G(t,t +T) exp( {7 b (v) dv) — 1 
exp( | bC) de) 
t,2((t — 7 (t))) | a(t)G(t,t H 
+Q (t, x(( (t))) | eG, t) Ga) 


1 
Lat rt a(t)G(t, t) — T 
+Q (t, x(( on( as) mo) 


f i ET s)) Q(s,x((s — 7 (s))) ds 
= Q(t, (t = T (4) 


(3.16) 
où E est donnée par (3.14). Enfin la substitution de (3.16) dans (3.15), achève 


la preuve. m 
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Lemme 3.5 Soit 


he [ow ii, BAT? emp Gf ma idi 


NW, 


Si 
A? > 4B, (3.17) 


alors, on a 


1 
2 
max{ f alu)du f vtt} < ; 
i 


Preuve. Comme a (t) + b(t) = p(t), b' (t) 


f odus f vidus f pijan et f oou f Wa 


Grâce à l'inégalité 


[reste [rome 


f or f Pore 


b 
où p(x) > 0, f p(x) dx > 0, inf f (x) > 0, nous arrivons à 


= exp 


f oa s00 (Ff mocan) , af a 200 (Ef nt 


Donc, nous obtenons aisément 


fe (u) af (u) du = [> (u) iu fa 


VV 
J 
bo 
io) 
cal 
gs 
TN 
le 
= 
ia 
5 
Le] 
O 
D 
a 
No 


T T T 
Sachant que | a (u) du + I b (u)du = f p(u) du, la conclusion découle 
0 0 0 


immédiatement de ce qui précède. m 
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Corollaire 3.2 Les fonctions G et E vérifient 


rs exp ( f pt) du) 


(1) | 


(3.18) 


3.3 Existence de solutions périodiques posi- 
tives 


Pour appliquer le théorème de point fixe de Krasnoselskii, nous avons besoin 
de définir un espace de Banach (B, ||.||), un sous ensemble convexe fermé D 
de B et de construire deux applications, l’une est une contraction et l’autre 
est continue et compacte. 

Soit (B, ||...) = (Pr, ||.) et D = {ye B:K <y< TL}, où K est une 
constante non négative et L est une constante positive. Nous exprimons 


l'équation (3.13) comme suit 


p (t) = (By) (t) + (Ay) (t) := (Hy) (4), 
où A, B : D — B sont deux applications définies par 
(Ay) (t) = f G(t, s) [F (s, hı (4 (8)) g1 (e(s—7(s)))) (3-19) 
+ (a (8) + a?(8)) Q (s, (ls = 7 (8)))] ds, 


et 


BAH = Qt pt- 70) (3.20) 
7 | p(s) E (t, s) Q (8, o((s— 7 (8))) ds. 
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Pour simplifier les notations, on introduit les constantes suivantes : 


rrexo( | p(w) di) 


e™ T 
Q = 5 = = ; 
(el — 1} =| Y (em — 1) (3.21) 
max | (é)|, = Ir(t)|, ma lp(t)| 


On suppose que la fonction Q(t, x) est localement lipschitzienne et continue 
par rapport à la variable x, donc il existe une constante strictement positive 


k telle que : 
Q(t, x) ~ Q(t, y)| < k læ E y|| ; vt € [0,7] , Vx E€ D. (3.22) 


Pour établir l’existence des solutions périodiques positives de (3.1), on dis- 


tingue deux cas Q(t, x) > 0 et Q(t,x) < 0 pour tout t E€ Ret x € D. 
Commençons d’abord par le premier cas où Q(t, x) > 0. 
On suppose qu’il existe une constante non négative kı et deux constantes 


positives kə et ø telles que : 


E (t,s) >c, V(t,s) € (0,7 | x [0, T], (3.23) 
kıx < Q(t, £) < kax, Vt € [0, T], Yz € D, (3.24) 

et 
ke <1, (3.25) 


ainsi, pour tout s € [0,7], z, y € D 
(1—k1) K + ABRTL 
yT 
< f(s h (2),g (9)) + (a(s) +0%(s)) Q(sy) (3.26) 
(1 = k2) L+ uokı T K 
< aT : 
Lemme 3.6 Supposons que les conditions (3.4)- (3-4), (3.7), et (3.23)- 


sont satisfaites, alors l'opérateur A défini par est compact. 


A 
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Preuve. 


Ligne directrice : En suivant les mêmes démarches que celles 
utilisées dans l’exemple 1.6 et en vertu du théorème d’Ascoli — Arzila, 


on peut montrer la compacité de l’opérateur À. 


Soit À l'opérateur donné par (3.19). Il est clair que A est continu, ainsi 
on peut montrer aisément que (Ay) (t + T) = (Ay) (t). 
Pour t € [0,T] et y, Y € D, ona 


ADOL < || GES) l6) a (l-7) 
+ (a (5) + a?(8)) Q (s, (ls = 7 (s)))] ds 
< ar k2)L + poh TK 


aT 
D’ot 
[Apl < (1 = k2) L + po TK, 
et par conséquent, l’ensemble A (D) est uniformément borné. 
Montrons maintenant que A (D) est équicontinu. Pour cela, soit y„ € D, avec 
n un entier naturel. On va calculer — (Ay,,) (t), puis on va montrer que cette 
dt 
dérivée est uniformément bornée. À l’aide de (3.2), (3.3) et (3.4), si on prend 
la dérivée par rapport à t dans (3.19), on obtient 
exp( f? a (v) dv) 
exp T a (v) dv) — 1) 
x [F (s, hi (Pn (8)) 91 (Pn (8 = T(s)))) 
—a (s) Q (5, 9,((s — T (s)))]. 


En utilisant (3.21) et (3.26), on arrive à 


d (Ag,) o {pe p) ARE ARMES 2p, 
dt aT 


SAP) = -oat s) + 
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pour certaine constante positive D. Alors la suite (Ay,,) est équicontinue. 
Le théorème d’Ascoli-Arzela assure l'existence d’une sous suite (Ay,,,) de 
(Apn) qui converge vers une fonction T—périodique et continue. Par consé- 
quent À est continu et A (D) est contenu dans un sous-ensemble compact de 


B. m 
Lemme 3.7 Supposons que la condition est satisfaite. Si 
(1+A8T)k <1, (3.27) 


alors l'opérateur B défini par est une contraction. 


Preuve. Soit B l'opérateur donné par (3.20). Il est clair que By est continu 


et il est aisé de voir que 
(By) (t +T) = (By) (+): 
En outre, pour tout t € [0,7] et y, y € D, ona 


(By) (t) — (By) (| 


< (Qt, P(t —7(t)))- Q (t, Y (t = 7 (¢)))| 
+f p(s) E (t, 5) [Q (s, 9 (s = T (s))) — Q (s, Y (s = 7 (s)))| ds 
< (AST) kle = tlla > 


il s'ensuit alors que 
Be — Byll < A+ ABT) k lly — tllo- 


Nous concluions par la condition (3.27) que B est une contraction. m 
L’équation intégrale obtenue (3.13) est une somme de deux applications ; 
une contraction et une application compacte et continue, donc dans ce cas 


on va appliquer le théorème de point fixe de Krasnoselskii. 
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Théorème 3.1 Supposons que les conditions (3.2)-(8.3), (8.7),[3.17 et (8.2)- 


sont satisfaites, alors l’équation admet au moins une solution 


positive T—périodique dans l’ensemble D. 


Preuve. 


Ligne directrice : La compacité de A et la contractivité de B sont 
assurées par les Lemmes|3.6 et|3.7, donc pour appliquer le théorème 


de Krasnoselskii, il nous reste qu’à vérifier sa première condition. 


D’après le Lemme A: D — B est compact et continu et d’après le 
Lemme B : D — B est une contraction. De plus, si y, Y € D, nous 
pouvons voir que 
(Ay) (t) + (BY) (t) 
t4+T 
= Qtut-ry)- | POEH- 7) 
t 


+f GESIE Meale- (9) 
+ (a'(s) + a(8)) Q (s, lls = 7 (8)))| ds 


t +T 


< by] Q (sols = 7 (5))) ds 
+a f TO C0) 
+ (a'(s) +%(s)) Q(s,9((s — 7 (s)))] ds 
< keb —pokiTK 4 je) “+ ee L, 
et 
(Ap) (D + (Be) (D 
> 


REA) @ (s, p(s = 7 (s))) ds 


hf [f (s, hi (9 (8), (Y (s —T(s)))) — a(s)Q (s, (8 — 7 (s)))] ds, 
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par conséquent 


(1 — ky) K + ABkoT L 
yT 


(Ay) (t) + (BY) (t) 


IV 


kK — XBRTL +7 
> K 


Maintenant, toutes les hypothèses du Théorème de Krasnoselskii sont 


satisfaites, alors il existe un point fixe x € D de l’opérateur 
H = A+B, 


l.e., 
x = Ar + Br. 
Enfin, le Lemme assure que ce point fixe est une solution de (3.1). m 


Dans le cas où Q(t, £) < 0, on suppose qu’il existe une constante négative ks 


et deux constantes non-positives k4 et o telles que : 
ksx < Q(t, x) < kaz, Vt € [0, T], Yz € D, (3.28) 


— k3 APT < 1, (3.29) 
ainsi, pour tout s € [0,T], z, y € D 
K (1 + kauoT) — kaL 
yT 
f (shi (2). (y)) + (a'(8) +) Q(s,y) (8:30) 
(1+ k3A8T) L — kak 
< aT . 
Théorème 3.2 Supposons que les conditions (8.2)-(8.3), (8.2), 
et (3.27)-(3.30) sont vérifiées, alors l’équation admet au moins une 


solution positive T— périodique dans l’ensemble D. 


IA 


Preuve. La preuve est analogue à la preuve du Théorème = 
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3.4 Récapitulatif 


Dans ce chapitre, nous avons abordé la question de l’existence des solu- 
tions périodiques positives d’une classe d'équations différentielles non li- 
néaires de type neutre et nous avons obtenu des résultats d'existence pouvant 
s'étendre aux EDN d’ordre supérieur à deux. Avant d’appliquer le théorème 
de point fixe de Krasnoselskii pour prouver l'existence des solutions, nous 
avons converti l'équation différentielle de type neutre en une équation in- 
tégrale s’écrivant comme somme de deux applications; la première est une 
contraction et l’autre est continue et est prouvée compacte moyennant le 


théorème d’Ascoli- Arzela. 
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CHAPITRE — 


Existence et unicité des solutions périodiques d’une 


classe d’équations différentielles itératives du second 


ordre 
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D” le but de prouver l’existence de solutions périodiques d’une classe 
d'équations différentielles du second ordre avec des termes itératifs, on 
a utilisé le théorème de point fixe de Schauder. Et via celui de Banach et 


sous certaines conditions appropriées, on a établi l’existence et l’unicité. 
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4.1 Introduction 


Les équations différentielles à retard dépendant de l’état sont utilisées 
pour modéliser mathématiquement plusieurs phénomènes dans différentes 
disciplines scientifiques comme l’électrodynamique, l’automatique, la biolo- 
gie, l’épidémiologie et l’économie, etc. L'étude de ces équations a commencé 
un peut avant avec le problème proéminent de deux corps en électrodyna- 
mique et qui a resté dans les limbes jusqu’à apparition des travaux de Driver 
[25] [28] 29] [80] dans les années 60 du siècle dernier. 

L’équation différentielle fonctionnelle à retard itérative qui est une rela- 
tion entre une fonction, ses dérivées et ses itérations, est un type particulier 
des équations différentielles à retard dépendant de l’état. 

Il semble que les travaux sur les équations avec termes itératifs ont débuté 
par un exemple simple mais très important traité par Babage en 1815 
où il est intéressé à trouver une fonction telle que sa n°" itération égale à 
la fonction elle même. 

Dans le siècle suivant, ce type d'équations a attiré l’attention de plusieurs 
chercheurs de l’époque et une très vaste littérature existe concernant l’étude 
de ce type d'équations par des techniques différentes. À titre d'exemples, on 
peut citer 


En 1968, P. Andrzej [5] a utilisé la méthode des approximations succes- 


sives de Picard, pour étudier l’équation 


w(t) = f (t,x (t) , x (x (t))) , 


où 0 est un point de domaine extrème à gauche. 
En 1984, Par le principe de l'application contractante, Eder [32] a obtenu 


quelques résultats d’existence d’une solution monotone unique de l'équation 


x' (t) = z (2 (t)), 
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avec æ(to) = to (to € [—-1,1]). 


L’équation 


avec x(a) = a, où a est un point de domaine extrême d’un intervalle bien 
défini, a été étudiée en 1990 par Wang [87] et par M. Fetkan [35] en 1993, 
où z(0) = 0. 

Quatre ans plus tard, la méme équation a été étudiée par W. Ge et Y. 
Mo [37 sous ’hypothése (to) = £o dans un intervalle compact donné, où 
les points extrêmes de l’intervalle sont deux points adjacents non nuls de f. 
Ce problème a attiré par la suite l’intension de Berinde qui a appliqué 
la méthode des opérateurs non expansives dans son travail publié en 2010. 
En 1995, Stanék [83] a étudié les propriétés globales des solutions de l’équa- 


tion différentielle itérative suivante : 
x’ (t) = a(t) +2 (x(t) 


et il a obtenu des résultats sur existence de solutions strictement monotones. 


En 1997, Si, Li et Cheng [79] ont considéré l’équation : 
2’ (t) =a) (6), 


et ils ont établi des conditions suffisantes pour assurer l’existence de solutions 
analytiques. Après une année, Si et Wang ont prouvé l'existence des 


solutions régulières de l’équation : 
a! (t) = Aya (E) + Aar (H) +. + Anr (t) + F). 


Plusieurs chercheurs se sont penchés ces dernières années sur l’étude de ce 
type d'équations en utilisant la technique du point fixe. Dans [92], H. Y. 
Zhaoa et J. Liu ont appliqué les théorémes de point fixe de Banach et de 
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Krasnoselskii pour établir l'existence, l’unicité et la stabilité de solutions 


périodiques de l’équation différentielle itérative du premier ordre suivante : 
Gull) = altt) + eH l (6) + + (ball (6) + F(t). 


Dans [91], H. Y. Zhaoa et M. Feckan ont utilisé le théorème de point fixe de 
Schauder pour montrer l'existence et la stabilité de solutions de l’équation 


suivante : 


z(t) = DCS Com (x! (0) + G0). 


m=14{=1 


Il faut signaler ici que la même dernière équation a été étudié par Liu et Si 
[56] en 2007 par la transformation de Schréder. 

Les travaux réalisés par la plupart des chercheurs jusqu’à maintenant res- 
tent concentrés sur les équations du premier ordre, par contre on ne trouve 
que très peu de travaux consacrés à l’étude des équations du second ordre. 
Dans ce chapitre et via les théorèmes de point fixe de Schauder et de Ba- 
nach, notre but majeur est d'établir des conditions suffisantes pour prouver 
l'existence et l’unicité des solutions périodiques de l’équation différentielle 
itérative du second ordre suivante : 


pat) + P(t) TOETO = “9 (¢, x(t), o7l(2),...,¢71(2)) 


+h (t, e(t), c?l(2),..., 24) (4.1) 


où p et q sont des fonctions positives continues, 7 = t, c(t) = z(t), 


sA) = x(a), s(t) = rP! (x(t)) et h, g : R! — R sont des 


fonctions continues par rapport à leurs arguments. 
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4.2 Préliminaires 


Pour tout T > 0 et L, M > 0, 


Pr(L,M) = {z € Pr, (xl < L, lx(t)-2(h) < M ltz- til, VibeR}, 


est un sous ensemble fermé, borné et convexe de Pr. Les fonctions h(t, £1, £2, … 
et g(t, £1, £2, ---, Zn) sont supposées T—périodiques en t et globalement lip- 


schitziennes en z1, ..., £n 1.6, 


h(t +T, 2£1,..., 8n) = h ie), g +T, 21,- En) = g (t, Z1, 3 En); 
(4.2) 
et qu’il existe n constantes positives k1, kə, ..., kn et n constantes positives 


C1, C2, .…, Cn telles que 


iti D ot oil 
i=1 
et | 
lg (t, Tı, yin) =g (t, Y1, …, ÿa)| < DC EZ = yill $ (4.4) 
i=1 


Lemme 4.1 [92/ Pour tout p, € Pr(L,M), 


o — pr 


n—1 
j=0 


Lemme 4.2 Supposons que les hypothèses (3.4), (8-H, et soient 


vérifiées. Si x € Pr(L, M), alors x est une solution de si et seulement 


si 


_ exp( f; b (v) dv) 
exp( f b(v) dv) — 1 


(4.6) 


83 


te) 


Chapitre 4. Existence et unicité des solutions périodiques d’une classe 
d'équations différentielles itératives du second ordre 


Preuve. 
Ligne directrice : En utilisant le Lemme 3.3 et certaines propriétés 


utiles de la fonction de Green citées dans le corollaire 3.1, on peut 


convertir l'équation (4.1) en une équation intégrale. 


Supposons que x € Pr (L, M) est une solution de (4.1). Du Lemme [3.3] dans 


le chapitre précédent, il est facile de voir que 


Comme 
[G(t, s)g (s, x(s), rPl(s), oe rl"l(s))|i = (), 


Il découle du Corollaire [3.1] que 
t+T d 
f G(t, s) a (s, 2(s), a? )(s), ..., 2l"!(s)) ds 
t 


= f g (s,æ(s),xP(s), …, al (s)) [E (t, s) — a(s)G(t, s)| ds. 


Il s’ensuit que 
6) = f [LE (t, s) — a(s)G(t, s)] g (s, x(s), cI(s), SI (s)) ds 


+ [ "alt s)h (t, x(t), tPl(t),..., oll (t)) ds. 


Ceci achève la preuve du Lemme. m 


paai) ds 
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4.3 Existence de solutions périodiques 


Dans cette section, on va utiliser le théoréme de point fixe de Schau- 
der pour prouver l'existence de solutions périodiques de l’équation (4.1). Du 


Lemme [1.2] on définit un opérateur S : Pr (L, M) — Pr par : 
t+T 
(Sy) (t) = 1 [E (t, 8) — a(s)G, 8)] 9 (sola) pPl(s), ..., p(s) ds 
t 
t+T 
+ f lE, 89h (5,019), 916), = 918) ds, (4.7) 
t 


et on va montrer qu'il est continu et compact et que Sy € Pr(L,M) pour 


tout y € Pr(L, M). Nous introduisons les notations suivantes : 


_ Texp (fp u) du) e™ = 
an D Ps V = exp (- b(v)do) 


1 
T M = max |a (t)| 
[expl a (u) du) — 1 lexp( fyb (u) du) — 1 telo,T] 
A2 = ae bl, a = ane ts 0,0,...,0)|, & = mig x 19 (6, 0,0, ..., 0)| 
€ 
n j=i-1 7 j=i-1 


G= tI) ok SoM, G=—,4+LS a Mi (4.8) 
i=1 j=0 i=1 j=0 
Lemme 4.3 Supposons que les conditions (3.9), (3.9), (3.7), et (4.4)- (4-4) 


sont satisfaites. Alors l’opérateur S : Pr (L, M) — Pr donné par est 


continu et compact. 


Preuve. 


Ligne directrice : La preuve est basée sur l’utilisation du 
théorème d’Ascoli — Arzila afin de montrer la compacité 
de Pr(L,M) et partant la compacité de S découle 


immédiatement de sa continuité. 
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Comme Pr (L, M) est un sous ensemble uniformément borné et équicontinu 
de l’espace des fonctions continues sur le compact [0,7], alors le théorème 
d’Ascoli— Arzila assure que Py (L, M) est un sous ensemble compact de cet 
espace. En outre, puisque tout opérateur continu envoie tout ensemble com- 
pact en un ensemble compact, il nous suffit de montrer que S est continu. 


Pour tout y, 0 € Pr(L,M), ona 


(Se) t) — ($8) (| 


$ [ lat) 1G, 8) |g (5, 806), 6%), 85) — g (s,6(5),6°%s),....6%(s)) ds 


+ FA 7 IG(t, s)| ln (s, p(s), pll(s), ..., p(s) — h (¢,0(s),0"(s), O) | = 


Sous les conditions (4.3), (4.4) et d’après le Corollaire [3.2] il vient 


Ip) (t) — (58) (| < BATY c p — 66 407 Y k: le" — g|), 
i=1 i=1 
En vertu du Lemme [4.1] il s’ensuit que 
n j=i-l 
(Se) (t) — (50) ©] <T X (B+ ad) ci + akı) D M? \lp— êl], 
i=l j=0 


ce qui prouve que l’opérateur S est continu. Donc, S est compact. m 


Lemme 4.4 Pour tout ti,t2€R 


ti+T 
| |G (62, 8) — G (t1, 8)|ds < |tz — tı| Tec [Tr (27 +1) +e" +1]. 
ti 
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Preuve. Soit tı,tə € R, on a 


ee — G (tı, s)| 
Longs eta oe en ftou+ are 
a aa du) — 1 ce b(u) du) — 
pet? exp fb) Face ae (v) do| du — [#7 exp feb ( v)du+ fit" a( (v) dv} du| 


lexp( (fo a u) du) — 1 [exol b (u) du) — 1 


f exp f b(vydo+ f a(o) d| du | exp | b(v)do+ f a(o) d| du 
t2 to u ty ty u 
f (elt COLE À ce du= i (el: bodo Le JE Hé) du 

t t 


2 
t2 i x 
-f (e b(v)dv of aoe) d 
t 


3 u u t2 u r 
_ i. efi a) v)dv i Siz OC b(v)dv efi k dr | (el: b(v)dv x elu aie) du 
t 


1 


t 
= ix JÈ a(v)dv of bv e (e AOLE eae) au- fo (cf b(v)dv eli atoa) dü 
t 


1 
S t 
| eliade y of bodo l ? Subd 6 fS.alv)de qn, 
t2 ty 
T t 
< eli b(u)do ela b(v)dv | fee a(v)dv x alo b(v)du du+ g left b(v)du x alo a(v)dv 
0 ty 
š F T t2 
< efa tde gle | edu f e"du 
0 tı 


< Te?" fef b(u)do gla b(v)dv 


< efi Odo _ Se, OC dv ee 


du 


=F e™ [to = ti] : 


Comme 


tı+T m s t2 
f fef b(u)du gla b(v)dv i alts b(v)dv 
ti 


ds 


tı +T : 
as -f eli b(v)dv 
tı 
T 
stiler (- f 17 iv) Ho tI, 
0 
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alors 


[ [x [eow Patan du- [exp [foarte Patan du 


T 
< Te?” ||b|| exp (- b(v) do) lto — t1| + Te? |t — til 
0 


= |tz — t| Te” (r \|b|| exp (- [ow do) + 1) (4.9) 


D’une manière similaire, on aura 


ds 


t2+T ti+T 
f ? als bod ofa T a(v)dv dy, _ f - alti LOUE aus a(v)dv day 


ti+T ti+T 
T ale b(u)dv o fa T aodo dy — I : eli b(v)de fiT a(v)de dy 


totT ù sir 
t1+T 


ti+T 447 
| | ale b(v)do free a(v)dv gay _ f | el blodo ofa T “ed 


< 


FT sir 
t 


1+T 


CEE p p 
= f (cls b(v)do ale nue edu a(v)dv dy, a 


ter 
= i 1 (cls b(v)dv ela ae Bigs alodo of Kodo gy 


UT u P 
Î efa MM ela” at qu 
t+T 


to+T i 
4 | . alta b(u)do ofa T a(v)dv dy, 
t 


1+T 


= alt b(u)du ela b(v)dv 


t1 +T 
f hia a(v)du of; b(v)dv dy, + 


titl 2T 2T 
| lef a(v)dv,, o b(v)dv 
s 


a T 2T 
+f lef b(u)dv, o a(v)du 
ti+T 


< fef b(v)dv alta b(v)du 


t2+T ù 
+r 
[7 aoten omy 
tı +T 


< efi b(v)dv = efa b(v)dv dü 


du 


Te + em [to = ti] : 
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Nous avons donc 


HT | plait . tere s 
‘a f ali oleo tied, f eft * vdv ofa +T dv du 
ty s s 
T 
< 27°e*" ||b|| exp (- [ b(v) do) ati] + Te?” |tz — til 
0 


T 
= |tz — t| T “a \|b|| exp (- b(v) iv) + on) (4.10) 
0 


Il s’ensuit de (4.9) et ( (4.10) que 
t1+T 
fO Cts) -Gl s)lds 


tı 


ds 


P i , IT 4m ||| exp (- fob -+ Te3m be e?™T al exp (- fro (v) do) X em 
> [42 ~ 41 
e E u) du) — 1 Je b (u) du) — 1 
T lloll exp (- h wen Dem + et) 4 ein + e2m 
= lt = t| 


lexp( Pa u) du) — 1 lexp( fyb (u) du) — 1 


< [te — ti| Te?™s [Thor we Hire Ls 


La preuve est achevée. m 


Lemme 4.5 Supposons que les conditions B-J- sont vérifiées. Si 


T(B+ahk)t+Taot<L, (4.11) 


et 
[ (2a + Tes (Trav (2e°" +1) +e” +. 1)) (A162 + Gy) + (28 + Trav) Go] < M, 


(4.12) 
alors 


S (Pr (L,M)) C Pp (L,M). 
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Preuve. 


Ligne directrice : Pour assurer l’inclusion ci-dessus, il suffit donc 
de montrer ces trois conditions : 

e S(p) € Pr, Yy € Pr(L, M) 

e |(Sp) (t2) — (S0) (te)| < M |t2 — til, Vti, te ER, Yy, 0 € Pr (L, M) 
e |S (p)|| < L, Ve € Pr (L, M). 


Il est aisé de voir que (Sy) (t +T) = (Sy) (t). Grâce au corollaire on 


arrive à 
t+T 
KSA O1 < (8+ a) f la (80C) PPE), = vli(s))] ds 
t 
t+T 
+ a | |h (t, p(t), PIG), ..., pI) | ds. 
t 
Sous les conditions (4.2)- (4.4) et en vertu du Lemme [4.1] on obtient 
Jr (s, 9(s), p(s), ..., o] 


= |h (s; 0(s), 92(s), .., pll(s)) — h(s,0,0,...,0)| 
+ |h (5,0,0, ...,0)| 


n  j=i-1 n  j=i-1 


<+) k D Mie <a +l k D MI < G, 
i=1 j=0 i=1 j=0 


et 
|g (s, p(s), pll(s), ..., o™(s)) | 
= |g (s, 0(s), p(s), ..., p™(s)) — g (s, 0,0, ..., 0)| 
+ |g (s,0,0,..., 0)| 
n ie n fei 
<po+ doc > M” || <p+LÝ G D M’ < Qo. 
i=1 j=0 i=1 j=0 
Donc 


(Sy) Œ| < T (8 + adı) Cı + Tags. 
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De (4.11), il vient 


[Syl] < L. 


Soit t1, t2 € R avec tı < to, ona 


(Sp) (t2) — (Sy) (ta) 


f Bedae 0(s),92l(s), --.9(s)) ds 


t2 


< 


= f E(t, 8) 9 (8, ¢(8), pPI (8), ~. wll(s)) ds 


+ f a(s)G (tz, s)g (s, p(s), pP! (s), ..., pll(s)) ds 


t2 


- 1  a(s)G(tr, 8)g (8, (8), PPIs), .., o"l(s)) ds 


ty 


+ f G(t2,s)h (s, p(s), p(s), ..., p(s) ds 


t2 


-f7 G(ti, s)h (s, p(s), p” (8), ..., pPl(s)) ds}. 


D’une part 


| ' E (ta, 8) g (5, (8), pPl(s),..., yl"l(s)) ds (4.13) 


t2 


- f Blt 8) 9(s, (8), es), s p8) ds 


< l “IE (ta, 8)| [9 (eple), ple), o"(s)) | as 


+ | Ets) |g (8, (8), 92s), … o"(s)) | ds 


1+T 


tı +T 
+f |E (tə, s) =E (t1, s)| |g (s, (5), pP (s), sapra) ds 
tı 
< 2 |t2 — ti] Co + [te — ti] TAgvc, 


= |t2 — t1] (28 + Thor) Ca, 
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et d’autre part 


tot T 
f a(s)G (t2, s)g (s, gls), pPl(s), ..., yl"l(s)) ds 


t2 


= f° als Gtts,8)g (8,000) ¢%3),.¢™)) ds 


tı 


< 


f daCi (8, 9(8), s) ds 


t2 


/  a(s)G(te, 8)9 (s, (s), pl(s), .., o"l(s)) ds 


f "abed — A(t, 8)] 9 (8, (8), e Od 


Du Lemme [4.4] et les conditions et (4.4), on aura 


f  E (ta, 8) g (8, (8), 9l(s), p! (8)) ds 


t2 


= f E (txs) g (8, 065), ets). p18) ds 


ti 


< Qadg [to — t1| Ca + Ar [te — t1| GT es (Thor (2er + 1) +e” + 1) 


= [ta = ti] Aio (2a + Lette (TAgv (2er + 1) + e™ + 1)) 5 


et 


92 


Chapitre 4. Existence et unicité des solutions périodiques d’une classe 
d'équations différentielles itératives du second ordre 


Grace au Lemme [4.4] et les conditions (2.5) —(1:6), on arrive à 


f G(t2,s)h (s, p(s), p"l(s), ..., wll(s)) ds 


t2 


-f G(t1, s)h (s, p(s), p(s), ..., pll(s)) ds 


ti 


< [lo — ti] ¢, (2a + Tes (Tr (2e°" +1) +e" +1)). (4.15) 


Ainsi, il s’ensuit de (4.13), (4.14) et (4.15) que 


(Sp) (t2) — (Sy) (t1) 
< [lo — ty] ((2a + Te (Thor (267 + 1) + e™ + 1)) (Arba + C1) + (28 + Thor) C3) . 


Moyennant la condition (4.12), on obtient 
(Sp) (t2) — (Se) (t1)| < M |t2 — ti]. 


Ce qui achève la preuve. m 


Voici maintenant notre résultat principal : 


Théoréme 4.1 Supposons que (3.4), (3.3), (3.7), (3.17), (4.4)- (4 
et sont vérifiées, alors l’équation admet au moins une solution 
x € Pr(L, M). 


Preuve. 
Ligne directrice : Le résultat découle immédiatement 


des Lemmes et 


Du lemme on voit que l’équation (4.1) admet au moins une solution 
x dans Pr (L, M) si et seulement si l’opérateur S défini par admet au 
moins un point fixe. Pour conclure, il suffit d’appliquer les Lemmes|4.3]et [4.5] 
ce qui prouve que les hypothéses du théoréme de Schauder sont satisfaites. 
Par conséquent, il existe au moins un point fixe x € Pr (L, M) de l’opérateur 


S. D'où l'existence d’au moins une solution de l’équation (4.1). m 
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d'équations différentielles itératives du second ordre 


4.4 Existence et unicité des solutions pério- 
diques 


Nous donnons pour clore ce chapitre un théorème qui assure l’existence 


et l’unicité des solutions de l’équation (4.1). 


Théorème 4.2 Supposons (8.2), (3.3), (8.1%), (4.9)- (4-4) sont satisfaites. 


Si 


j=i-1 


ry B+ah)c+ak;) D. Mi <1, (4.16) 


j=0 


alors l’équation admet une solution unique x € Pr(L, M). 


Preuve. Soit y, 0 € Pr(L, M). Par une démonstration analogue à celle du 


Lemme on obtient 


(Sy) (t) — (58) (t Box Gao) So M’) heal 


i=1 j=0 

ce qui entraîne que moyennant l’hypothèse (4-16), S est une contraction. 
Sous la condition citée ci dessus, et en vertu du théorème de point fixe de Ba- 
nach, S admet un point fixe unique dans Pr (L, M), ce qui prouve l'existence 


d’une solution unique de l’équation (4.1). m 
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4.5 Récapitulatif 


Nous avons abordé dans ce chapitre quelques questions liées à l’existence, 
l’unicité et la périodicité des solutions d’une classe d'équations différentielles 
fonctionnelles itératives du second ordre. Les résultats sont obtenus en se 
basant essentiellement sur une méthode combinant le théorème de point fixe 
de Schauder et les propriétés de la même fonction de Green du chapitre 
précédent. 

Nous montrons aussi que, sous certaines conditions et en vertu du théo- 
rème de point fixe de Banach, la même équation différentielle itérative pos- 


sède une solution unique. 
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| Conclusion générale et Perspectives 


L’intérét de la théorie de point fixe pour l’étude des équations fonctionnelles 
non linéaires a attiré mon attention et c’est pour cette raison que je l’ai choisi 
comme sujet de ma thèse de Doctorat où nous sommes intéressés à établir 
des conditions suffisantes pour assurer l’existence, l’unicité, la périodicité et 
la positivité des solutions de deux classes d'équations différentielles à retard 
variable du second ordre où la première classe, est une classe d'équations de 
type neutre tandis que la deuxième est d'équations de type itératif. 

Au cours de ce travail, nous avons présenté un état de l’art sur les équa- 
tions fonctionnelles à retard et nous avons également mis en évidence la 
grande efficacité de la technique qui combine quelques propriétés utiles des 
fonctions de Green de certaines EDO linéaires avec les théorèmes de point 
fixe. 

Cette technique puissante repose sur l’inversion de l’équation considérée 
en une équation intégrale avant d’appliquer le théorème de point fixe de Kras- 
noselskii ou celui de Schauder combiné avec quelques propriétés d’une fonc- 
tion de Green d’une équation différentielle linéaire du second ordre. En outre, 
nous avons appliqué le principe de l’application contractante de Banach sur 


l'équation intégrale inversée et sous certaines contraintes déterminées, on a 
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Conclusion et perspectives 


démontré l’existence et l’unicité de la solution périodique de l’équation dif- 
férentielle de type itératif. 

Ce type d'équations est à la mode aujourd’hui et un changement soit 
dans l’espace de départ ou soit le type du retard ouvre la voie vers pleusieurs 
autres pistes interessantes ayant trait au points suivants : 

— Il apparaît intéressant de démontrer la stabilité de ces solutions. 

— Il serait aussi très utile de pouvoir étendre l’étude aux équations diffé- 
rentielles à retard d’ordre supérieur à deux, équations différentielles à retard 
distribué ou équations fractionnaires à retard. 

— Il semble importent d'étudier l’existence de solutions presque pério- 


diques ou anti périodiques. 
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| Annexe 


I. Résolution des EDR par MATLAB 
Les principales fonctions de MATLAB qui permettent la résolution des équa- 


tions différentielles à retards sont dde23, ddesd et ddensd. 


1. Résolution d’une équation différentielle A retard constant par 
dde23 
Le solveur dde23 de Matlab, basé sur les méthodes de Runge-Kutta 
explicite est développé par S. Thompson et L. Shampine pour résoudre 
les équations différentielles avec des retards constants. 


Syntaxe de base 


sol dde23 ( ddefile, lags, history, tspan) 


sol dde23 ( ddefile, lags, history, tspan, options), 


ot les arguments d’entrée sont définis par : 

ddefile : fonction au second membre de PEDR. 

lags : vecteur des retards constants positives. 

history : fonction historique, c’est un vecteur colonne donnant les 
valeurs initiales des solutions. 


tspan : intervalle d’intégration (temps ot on veut calculer la solution). 
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2. Résolution d’une équation différentielle 4 retards variables dé- 
pendant continiment du temps ou de l’état par ddesd 
La fonction ddesd permet la résolution des équations différentielles avec 
des retards variables. 


Syntaxe de base 


sol = ddesd (ddefun, delays, history, tspan) 


sol ddesd (ddefun, delays, history, tspan, options), 


où les arguments d’entrée sont définis par : 

ddefun : fonction créée qui évalue le second membre de PEDR 

delays : vecteur colonne qui représente les retards qui peuvent dé- 

pendre du temps ou de l’état. 

history : fonction historique, c’est un vecteur colonne donnant les va- 

leurs initiales des solutions. 

tspan : intervalle d’intégration (temps où on veut calculer la solution). 
3. Résolution d’une équation différentielle à retard de type neutre 

par ddensd 


La fonction ddensd permet la résolution des équations différentielles de 
type neutre avec des retards indépendants de l’état. 


Syntaxe de base 


sol — ddensd (ddefun, dely, delyp, history, tspan) 


sol = ddensd (ddefun, dely, delyp, history, tspan, options). 


où les arguments d’entrée sont définis par : 
ddefun : fonction créée qui évalue le second membre de EDR 
dely : vecteur colonne qui représente les retards. 


delyp : vecteur colonne qui représente les retards dans les dérivées. 
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history : fonction historique, c’est un vecteur colonne donnant les 
valeurs initiales des solutions. 


tspan : intervalle d’intégration (temps ot on veut calculer la solution). 


Un tutoriel sur la résolution des équations différentielles à retard via 


MATLAB par Shampine et Thompson peut être trouvé sur le site 
Webhttp ://www.runet.edu/thompson/webddes/ 


II. Résolution des EDR par Maple 

La premiére procédure est appelée fhmpas, repose sur une résolution 
exacte des EDO successives. Elle devient lente dés que le nombre d’itérations 
devient grand et doit étre alors remplacée par des procédures de résolution 
approchée, appelées fhmpasn et fhmpasndp, utilisant les ressources de option 


numeric dans dsolve. 
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